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1.0 INTRODUCTION 


i - 


National Aeronautics and Space Administration (NASA) Contract NAS7-961 
[A Small Business Innovation and Research (SBIR) contract from NASA] 
involved research dealing with remote structural damage detection using the 
concept of substructures [1], Several approaches were developed; the main 
two were: (1) the module (substructure) transfer function matrix (MTFM) 
approach and (2) modal strain energy distribution method (MSEDM) . Either 
method can be used with a global structure; however, the focus was on sub- 
structures. 

As part of the research contract, computer software was to be developed 
which would implement the developed methods. This was done and it was used 
to process all the finite element generated numerical data for the research. 
The software was written for the IBM AT personal computer. Copies of it 
were placed on floppy disks (see Section 4.0). 

This report serves as a user's manual for the two sets of damage 
detection software. Sections 2.0 and 3.0 discuss the use of the MTFM and 
MSEDM software, respectively. 
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2.0 MODULE TRANSFER FUNCTION MATRIX METHOD 


This method/approach is discussed in detail in Section 3.0 of Reference 
1. The corresponding software consists of several programs. To discuss 
this, it is easiest to describe the process followed in doing the related 
part of the research. 

The basic steps followed are given as follows: 

a) obtain transient response of finite element model corresponding 
to undamaged and damaged states; this was done using COSMOS/M 
[ 2 ]; 

b) process transient response data to obtain corresponding data in 
MAC/RAN IV [3] format; 

c) run program for determining the module transfer function matrix 
[MTFM] for the undamaged and damaged states; and 

d) determine differences between the MTFM for the two states. 

The above steps were followed for each damage level (state) for a given 
model . 

The MTFM software was built around the computer code MAC/RAN IV. This 
code is a well developed digital time series analysis package. The most 
important part of it that is related to the subject research is its module 
(subroutine) TRANS; it is a multi-input, single-output (MISO) code. 
However, the research required a multi-input, multi-output (MIMO) code. To 
deal with this, a computer code was written that would execute TRANS for 
each required output (one at a time) — the code looped through TRANS the 
number of times needed. The TRANS output files were merged together to form 
a rectangular matrix corresponding to the MIMO problem. 

The details of how to run the MTFM and damage indicator computer codes 
are described in Table 2.1. All the codes referenced in this table can be 
found on each of the 5-1/4-in. floppy disks entitled, "Structural Module 
Integrity Software, ANCO Engineers, Inc., Report 1311.05, August 1988, NASA 
Contract NAS7-961." Table 2.2 describes the contents of the floppies that 
is related to the MTFM method. Appendix A contains a listing of the 
contents of a floppy, including the source listings. 

An example of an MTFM and damage indicator run is given in Table 2.3. 
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TABLE 2.1: RUNNING THE MTFM AND DAMAGE INDICATOR COMPUTER COOES 


i 


The codes were written to be used in connection with the finite element 
computer code COSMOS/M. The sequence followed in running the software 
package is given as follows: 



Description 


Start process using COSMOS/M output time histories corres- 
ponding to the undamaged and damaged substructure states. 


Use macro NASAMAC.BAT - this macro reads the COSMOS/M 
output time histories and outputs time histories in MAC/RAN 
IV format. 


To do: 


2-1 Create NASAMAC.BAT -- 

INTRFS.EXE 

MAC423 /R 10000 <FRMCHK.SCF >MAC0UT.0UT 

DEL MAC0UT.0UT 

NASADAM 

2-2 Create NASADAM — 


NASA.EXE * 
NASAD 


* NASA.EXE creates NASAD . BAT , NASAREN.BAT and MATRX.DAT. 

2-3 Run NASAMAC.BAT 

Input to INTRFS.EXE is given as follows: 

a) Enter COSMOS plot filename to process 

b) Enter output MACRAN.SCF filename - include extension 
• SCF 

c) Enter output filename MACRAN.SDF filename - include 
extension . SDF 

d) Enter the number of points in the time history 

e) Do you want to truncate the number of points? 

f) Enter the number of loads to use from COSMOS 

g) Is there load data in the COSMOS plot file? 

h) Enter the number of accelerations to use from COSMOS 
(this must be the total number of node points in the 
COSMOS output file) 

i) Enter channel description for each acceleration used, 
e.g., — >N0DE1,TX; — >N0DE6,TX; — >N0DE9,TX 
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TABLE 2.1 (continued) 


0 Step Description 

3 Input for NASA . EXE — 

This program is used to create a batch file which is used 

to execute MAC/RAN once for each output channel*. 

Input is given as follows: 

a) Enter input MACRAN FILENAME - including extension .SDF 

b) Number of points (time history points) 

c) Number of points/segments (the number of segments refers 
to how many segments the time histories are broken up 
into; for this study four (4) segments were used) 

d) Enter output filename - no file extensions please 

e) Enter the number of input channels (<=10) (this is the 
number of inputs to use for the calculation of the 
transfer function matrix) 

f) Specify MAC/RAN channel to use for each input (in 
MAC/RAN, the first channel is time (values of time in 
time history); the remaining channels correspond to the 
node numbers and their order in the COSMOS output file, 
e.g., time. Node 1, Node 6 , and Node 9 correspond to 
MAC/RAN Channels 1, 2, 3 and 4, respectively. 

g) Specify MAC/RAN channel for each output 

h) Is this a (d)amaged or (u)ndamaged run? 

(Enter "D” or "U") 

i) Enter a descriptive title for this run 


* MAC/RAN is a multi-input, single-output transfer function 
processor code. 

4 At this point in the process, a multi-input, multi-output 

transfer function matrix has been created for the substruc- 
ture. One matrix is created for each run of INTRFS.EXE, 
either undamaged or damaged condition. A matrix for both 
damaged states must be obtained before proceeding with the 
analysis. The elements of the MTFM are plotted by editing 
the MTFM data file to obtain the needed data for the 
element, and then by using the plotting features in MAC/RAN 


2-3 



i 


TABLE 2.1 (continued) 




The final step in this process is to calculate the damage 
indicators. The results are plotted using the plotting 
features in MAC/RAN 
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TABLE 2.2: CONTENTS OF FLOPPY DISK THAT ARE RELATED 

TO THE MTFM METHOO 
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TABLE 2.2 (continued) 
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TABLE 2.3: EXAMPLE OF MTFM METHOD 


¥ 


The input to this set of programs is rather lengthy, therefore, the input 
and output listings for the problem discussed in Subsection 4.2.2 of 
Reference 1 is shown below. The COSMOS/M time thistory plot files used are 
FRMWWUX.PLT and FRMWWD2X.PLT for the undamaged and damaged states, respect- 
ively. 


INTRFS - PROGRAM TO REFORM COSMOS FILE TO MACRAN FILE 
VIF'AC-VIBRATION PACKAGE REVISION NO.: 1.00 

ANCO ENGINEERS, INC. <C) Copyright 1904,1985,1986 

ENTER COSMOS PLOT FILENAME TO PROCESS : FRMWWUX.PLT 


ENTER 

OUTPUT 

MACRAN 

.SCF 

FILENAME - 

INCLUDE 

.SCF : FRMWWUX . SCF 

ENTER 

OUTPUT 

MACRAN 

. SDF 

FILENAME - 

INCLUDE 

.SDF : FRMWWUX. SDF 


ENTER THE NUMBER OF POINTS IN THE TIME HISTORY : 1024 

DO YOU WANT TO TRUNCATE THE NUMBER OF POINTS? N 

ENTER THE NUMBER OF LOADS TO USE FROM COSMOS : 0 

IS THERE LOAD DATA IN THE COSMOS PLOTFILE? : N** 


ENTER THE NUMBER OF ACCELERATIONS TO USE FROM COSMOS : 3 

ENTER CHANNEL DESCRIPTION FOR EACH ACCELERATION USED 
— > NODE .1 , TX 
— > N0DE6 , TX 
— > N0DE9 , TX 

FILE HAS A TIME STEP OF : 2.000E-02 SECONDS 

READING IN ACCELERATION DATA 
REFORMATTING TO COSMOS -.SCF FORMAT 


Stop - Program terminated. 

21: 14: 40 
E: \ 
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TABLE 2.3 (continued) 


21: 14*40 
E: \ 

>TYPE FRMWWUX . PLT 


1024 9 tim* d 1 

1 V 1 

1 a 1 1 d 1 6 v 1 

6 a 1 

6 d 1 

1 9 a 1 + 

. 20000E-01 . OOOO 

0000 .0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

. 40000E-0 1 .0000 

0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

. 60000E-01 . 0000 

0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

. 80000E-01 .0000 

0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

. 1 OOOOE+OO . 0000 
0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

. 1 2000 . 0000 
0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

.14000 .0000 

0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

. 1 6000 . 0000 
0000 . 00 00 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

. 1 8000 . 0000 
0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

. 20000 . 0000 

1 . 9800 . 0000 

0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

2. 0000 . 0000 
0000 . 0000 

. 0000 
. 0000 

. 0000 

. 0000 

. 0000 

. 0000 

2. 0200 , 0000 
3388E-04 3388E-02 

. 0000 

4687E-01 

.2919 

2034E-05 

2034E-03 

. 2716 

2- 0400 - 0000 

2285E-03 -.1607E-01 

. 0000 

-.5564E-03 

. 9292 

1372E-04 

965 1 E— 03 

. 8734 

2. 0600 . 0000 
7955E-03 — . 4063E-0 1 

. 0000 
. 3654 

1 . 891 

4776E-04 

2439E-02 

1 . 799 

2- 0800 . 0000 
1 850E-02 - . 6479E-0 1 

. 0000 
1 . 392 

2. 283 

-.11 1 IE-03 

3890E-02 

2. 230 

2- 1000 .0000 
2994E-02 -.4961E-01 

. 0000 
2.985 

. 5753 

1 798E-03 

2979E-02 

. 7200 

2 . 1 200 . 0000 
3197E-02 . 2930E-01 

. 0000 

T crcn 
•Jt . JJ*i. 

-1 . 930 

1919E-03 

. 1 759E-02 

- 1 . 600 

2.1400 .0000 

1676E-02 .1228 

. 0000 

1 . 621 

-2.245 

-. 1006E-03 

. 7372E-02 

-2 . 0 1 3 

2.1600 .0000 
1020E-02 .1468 

. 0000 

-1.779 

-.3218 

. 6125E-04 

. 8818E-02 

~. 4093 

2.1800 .0000 
3356E-02 .8670E-01 

. 0000 
-3.964 

. 5932 

. 2015E--03 

. 5206E-02 

. 3 1 96 

2. 2000 . 0000 

4280E -02 . 56B5E-02 

2. 2200 . 0000 

. 0000 

-3.991 

. 0000 

4468 

. 2570E-03 

. 3414E-03 

-.6596 
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TABLE 2.3 (continued) 


21s 16:42 
Es\ 

>TYPE FRMWWUX . 3CF 

JBHEAD TIME HISTORIES FROM COSMOS FILE - FRMWWUX. PLT 

FI LOUT FRMWWUX. SDF 

CTAPE NASA 1311.05 

CTAPE COSMOS/M 4. 

CTAPE TIME 
CTAPE NODE 1 , TX 
CTAPE N0DE6 , TX 
CTAPE N0DE9 , TX 


CTAPE 

. 000 E +00 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

2 . 0 COE— 02 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

4 . 000E-02 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

6 . 000E-02 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

8- 000E-02 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 . 000E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 - 200E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 - 400E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 . 600E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 . 800E-01 

. OOOE+OO 

. OOOE+OO 

.OOOE+OO 

CTAPE 

2. 000E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

2.200E-01 

. OOOE+OO 

. OOOE+OO 

.OOOE+OO 

CTAPE 

2. 400E-01 




CTAPE 

1 . 900E+00 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 . 920E+00 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 - 940E+00 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 . 960E+00 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 , 980E+00 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

2. OOOE+OO 

2. 919E-01 

2. 716E-01 

-4. 687E-02 

CTAPE 

2- 020E+00 

9.292E-01 

8.734E-01 

-5 . 564E-04 

CTAPE 

2. 040E+00 

1 . 891E+00 

1 . 799E+00 

3. 654E-01 

CTAPE 

2 . 060E+00 

2. 283E+00 

2. 230E+00 

1 . 392E+00 

CTAPE 

2. 080E+00 

5. 753E-01 

7. 200E-01 

2. 985E+00 

CTAPE 

2- 100E+00 

- 1 . 930E+00 

-1 . 600E+00 

3 . 552E+00 

CTAPE 

2. 120E+00 

-2. 245E+00 

-2. 013E+00 

1 . 621E+00 

CTAPE 

2, 140E+00 

—3. 21 BE— 01 

-4.093E-01 

-1 . 779E+00 

CTAPE 

2. 160E+00 

5. 932E-01 

3. 196E-01 

-3. 964E+00 

CTAPE 

2. 180E+00 

-4. 468E-01 

-6.596E-01 

-3. 991E+00 

CTAPE 

2. 200E+00 

-2. 506E-01 

-4. 247E-01 

— 3. 1 50E+OO 

CTAPE 

2, 220E+00 

1 . 846E+00 

1 . 631E+00 

-1 . 721E+00 

CTAPE 

2. 240E+00 

2. 696E+00 

2. 603E+00 

1.1 63E+00 

CTAPE 

2 . 260E+00 

1 . 396E+00 

1 . 594E+00 

4 . 682E+00 

CTAPE 

2, 2B0E+00 

-3. 275E-01 

9. 23 IE— 02 

6 . 664E+00 

CTAPE 

2 , 3 00 E +00 

- 1 . 079E+00 

-6. 677E-01 

5 . 774E+00 

CTAPE 

2. 320E+00 

-3. 21 IE— 01 

-1 . 588E-01 

2. 383E+00 

CTAPE 

2. 340E+00 

1 . 079E+00 

9. 207E-01 

- 1 . 565E+00 

CTAPE 

2. 360E+00 

1 . 284E+00 

9. 743E-01 

-3. 870E+00 

CTAPE 

2, 380E+00 

2. 571E-01 

1 
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TABLE 2.3 (continued) 


{ INTRFS - PROGRAM TO REFORM COSMOS FILE TO MACRAN FILE 
VIPAC-VIBRATION PACKAGE REVISION NO.: 1.00 

ANCO ENGINEERS, INC. (C) Copyright 1984,1985,1986 

ENTER COSMOS PLOT FILENAME TO PROCESS : FRMWWD2X . PLT 

ENTER OUTPUT MACRAN . SCF FILENAME - INCLUDE . SCF » FRMWWD2X . SCF 

ENTER OUTPUT MACRAN . SDF FILENAME - INCLUDE . SDF : FRMWWD2X . SDF 

ENTER THE NUMBER OF POINTS IN THE TIME HISTORY : 1024 

DO YOU WANT TO TRUNCATE THE NUMBER OF POINTS? N 

ENTER THE NUMBER OF LOADS TO USE FROM COSMOS : O 

IS THERE LOAD DATA IN THE COSMOS PLOTFILE? : N** 


ENTER THE NUMBER OF ACCELERATIONS TO USE FROM COSMOS : 3 

ENTER CHANNEL DESCRIPTION FOR EACH ACCELERATION USED 
— > N0DE1 , TX 
— > N0DE6 , TX 
— > N0DE9 , TX 

FILE HAS A TIME STEP OF : 2.000E-02 SECONDS 

READING IN ACCELERATION DATA 
REFORMATTING TO COSMOS -.SCF FORMAT 


Stop - Program terminated. 

21 : 27: 32 
E: \ 
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TABLE 2.3 (continued) 


, E:\ 

; TYPE FRMWWD2X . PLT 


1024 

9 time d 1 

1 V 1 

1 a 1 

d 1 6 v 1 

6 a 1 

6 d 1 

1 9 

a 1 9 






. 20000E 

-01 .0000 

. 0000 

. OOOO 

. 0000 

. OOOO 

. OOOO 

0000 

. 0000 

. 0000 





. 40000E 

-01 .0000 

. 0000 

. OOOO 

. oooo 

. oooo 

. oooo 

0000 

. 0000 

. 0000 





. 60000E 

— 0 1 . 0000 

. 0000 

. oooo 

. oooo 

. oooo 

. oooo 

0000 

. 0000 

, 0000 





. 90000E 

-01 .0000 

. 0000 

. oooo 

. 0000 

. oooo 

. oooo 

0000 

. 0000 

. 0000 





- 10000E+00 .0000 

■ 0000 

, oooo 

. OOOO 

. oooo 

. oooo 

0000 

. 0000 

• 0000 





. 12000 

. 0000 

. 0000 

. oooo 

. oooo 

. oooo 

. oooo 

0000 

, 0000 

. 0000 





. 14000 

- 0000 

. 0000 

. oooo 

. oooo 

. oooo 

. oooo 

0000 

. 0000 

- 0000 





. 1 6000 

. 0000 

. 0000 

. oooo 

. oooo 

. oooo 

. oooo 

0000 

, 0000 

. 0000 





. 1 8000 

. 0000 

. 0000 

. oooo 

. oooo 

. oooo 

. oooo 

0000 

• 0000 

. 0000 





. 20000 

. 0000 

. 0000 

. oooo 

. oooo 

. oooo 

. oooo 

0000 

. 2 

. 0000 

, 0000 





1 . 9800 

. 0000 

. 0000 

. oooo 

. oooo 

. oooo 

. oooo 

0000 

- 0000 

. 0000 





2. 0000 

- 0000 

. oooo 

. oooo 

. oooo 

. oooo 

. oooo 

0000 

. 0000 

. 0000 





2 . 0200 

. 0000 

. 0000 

. 2919 

-. 1 384E-05 

-. 1384E-03 

. 2781 

3239E-04 

-.3239E-02 

3205E-01 





2. 0400 

. 0000 

. 0000 

. 9292 

9448E-05 

6680E-03 

. 8901 

22 1 1 E— 03 

1563E-01 

. 1368E-01 





2. 0600 

. 0000 

. 0000 

1 . 891 

3343E-04 

-. 1 730E-02 

1 . 824 

7825E-03 

4050E-01 

. 3196 





2. 0800 

. 0000 

. 0000 

2.283 

7955E-04 

-. 288 IE— 02 

2.236 

1862E-02 

- . 6744E-0 1 

1 . 161 





2. 1000 

. 0000 

. 0000 

■ 7 <*j o 

-. 1 33QE-03 

-. 2546E-02 

. 6568 

3132E-02 

-.5959E-01 

2.483 





2. 1200 

. 0000 

. 0000 

- 1 . 930 

155 IE -03 

. 4194E-03 

-1.714 

3630E-02 

. 9818E-02 

3. 103 





2 „ 1400 

. 0000 

. 0000 

-2. 245 

-. 1076E-03 

. 4327E-02 

-2.069 

2519E-02 

. 1013 

1 .869 





2. 1600 

. 0000 

. 0000 

-.3218 

4676E-05 

. 5967E-02 

— • 3335 

1094E-03 

. 1397 

- . 5962 





2. 1800 

. 0000 

. 0000 

. 5932 

. 1005E-03 

. 4547E-02 

. 4629 

235 IE-02 

. 1064 

-2.457 





2- 2000 

. 0000 

. 0000 

4468 

. 1673E-03 

. 2136E-02 

5576 

3916E-02 

. 4999E-01 

-3. 040 





2 . 2200 

. 0000 

. oooo 

2506 

. 1 866E-03 

-. 2003E-03 

3734 
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TABLE 2.3 (continued) 


21:29:03 
E: \ 


>TYF'E 

FRMWWD2X. SCF 





JBHEAD 

TIME HISTORIES 

FROM COSMOS FILE 

- FRMWWD2X . 

PLT 


FI LOUT 

CT APE 

CTAPE 

CTAPE 

CTAPE 

CTAPE 

CTAPE 

CTAPE 

FRMWWD2X.8DF 
NASA 1311.05 
COSMOS/M 
TIME 

NODE 1 , TX 
N0DE6 , TX 
N0DE9 , TX 

4. 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


2. 000E-02 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


4. 000E-02 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


6 . 000E-02 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


8. 000 E- 02 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


1 . OOOE-Ol 

. OOOE+OO 

. OOOE+OO 

• OOOE+OO 

CTAPE 


1 . 200 E—0 1 

. OOOE+OO • 

. OOOE+OO 

. OOOE+OO 

CTAPE 


1 . 400E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


1 . 600E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


1 . 800E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


2 . OOOE-O 1 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 


2. 200E-01 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 


CTAPE 

1 . 920E+00 

• OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 . 940E+00 

- OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 . 960E+00 

. OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

1 . 980E+00 

• OOOE+OO 

. OOOE+OO 

. OOOE+OO 

CTAPE 

2. OOOE+OO 

2. 919E-01 

2. 78 IE- 01 

-3. 205E-02 

CTAPE 

2. 020E+00 

9.292E-01 

8.901E-01 

1 . 368E-02 

CTAPE 

2. 040E+00 

1 . 891E+00 

1 . 824E+00 

3. 196E-01 

CTAPE 

2. 060E+00 

2. 283E+00 

2 . 236E+00 

1 . 161E+00 

CTAPE 

2. 080E+00 

5. 753E-01 

6. 568E-01 

2. 483E+00 

CTAPE 

2. 100E+00 

-1 . 930E+00 

-1 . 714E+00 

3. 103E+00 

CTAPE 

2. 120E+00 

-2. 245E+00 

-2. 069E+00 

1 . 869E+00 

CTAPE 

2. 140E+00 

-3. 218E--01 

-3. 335E-01 

-5. 962E-01 

CTAPE 

2. 160E+00 

5. 932E-01 

4. 629E-0 1 

-2. 457E+00 

CTAPE 

2. 180E+00 

-4. 468E-01 

-5. 576E-01 

—3 , 040E+00 

CTAPE 

2. 200E+00 

-2. 506 E-0 1 

-3. 734E-01 

-3. 126E + 00 

CTAPE 

2 . 220E+00 

1 . 846E+00 

1. . 654E+00 

-2. 643E+00 

CTAPE 

2. 240E+00 

2. 696E+00 

2. 551E+00 

-6. 969E-01 

CTAPE 

2 . 260E+00 

.1 . 396E+00 

1 . 433E+00 

2 . 265E+00 

CTAPE 

2. 280E+00 

-3. 275E-01 

~! . 154E-01 

4. 637E+00 

CTAPE 

2 . 300E+00 

-1 . 079E+00 

-8. 1 04 E—0 1 

5. 213E+00 

CTAPE 

2. 320E+00 

-3.21 IE-01 

-1 . 4 28 E— 01 

3. 853E+00 

CTAPE 

2. 340E+00 

1 . 079E+00 

1 .099E+00 

1 . 537E+00 

CTAPE 

2. 360E+00 

1 . 284 E +00 

1.21 IE +00 

-4. 090E-01 

CTAPE 

2. 380E+00 

2. 571E-01 

1 „ 790E-01 

-1 . 57 IE +00 
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TABLE 2.3 (continued) 


NASA - TO GENERATE A MACRAN BATCH FILE 
VIPAC-VIBRATION PACKAGE REVISION NO.: 1.00 

ANCO ENGINEERS, INC. (C> Copyright 1984,1905,1986 


ENTER INPUT MACRAN FILENAME (INCLUDING . SE^) : FRMWWD2X . SDF 

NUMBER OF POINTS : 1024 

NUMBER OF POINTS/SEGMENT ; 256 

ENTER OUTPUT FILENAME - NO FILE EXTENSIONS PLEASE s FRMWD2 
ENTER THE NUMBER OF INPUT CHANNELS ( <= iO ) : 1 

SPECIFY MACRAN CHANNEL TO USE FOR EACH INPUT 
— > 3* 


SPECIFY MACRAN CHANNEL FOR EACH OUTPUT 

— > 4 

IS THIS A (D)AMAGED OR ( U ) NDAMAGED RUN? :D 

ENTER A DESCRIPTIVE TITLE FOR THIS RUN : FRAMWB DAMAGED, XBMO T 


Stop - Program terminated. 

2:24: 23 
: \ 
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3.0 M00AL STRAIN ENERGY DISTRIBUTION METHOD (MSEOM) 


This method (modal strain energy distribution method) is discussed in 
detail in Section 6.0 of Reference 1. The corresponding software consists 
of several programs. The basic steps followed in using it are given as 
fol lows: 

a) run a program to determine the strain energy in each finite 
element for each mode; 

b) determine the level of damage associated with each mode; this 
can be done by using the damage indicator discussed in 
Reference 1 ; and 

c) run program to determine possible locations of damage. 

The details of how to run the MSEDM computer codes are described in 
Table 3.1. All of the codes referenced in this table can be found on each 
of the 5-1/4-in. floppy disks entitled, "Structural Module Integrity 
Software, ANCO Engineers, Inc., Report 1311.05, August 1988, NASA Contract 
NAS7-961 . " Table 3.2 describes the contents of the floppies that are related 
to the MSEDM. Appendix A contains related source code listings. 


An example of an MSEDM run is given in Table 3.3. 



i 


TABLE 3.1: RUNNING THE MSEDM COMPUTER CODES 


The following three computer codes must be run in using the MSEDM. The 
required inputs are listed for each code. 


MODSTRN - This program computes element modal strain energies for 
a simply-supported beam due to bending deformation only. The user 
must input the number of elements and modes to use. The output 
element modal strain energies are contained in the file "MS.DD". 
This file is a binary file. 

DAMIO - This program allows the user to input damage parameters to 
be used in damage assessment. Damage parameters are input on a 
mode-by-mode basis. The user also specifies the output file con- 
taining the damage parameters. This output file is an ASCII file. 

MODSUM - This program performs the damage assessment based on the 
modal strain energies and the damage parameters. The user must 
input the file name containing the modal strain energies and the 
damage parameters. The user also specifies the damage fraction 
(g) used in assessing damage. Alternately, the user can select 
the option of summarizing the element strain energies for a given 
mode. In either case, output data is contained in the file 
"MODSUM. DD". 
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TABLE 3.2: CONTENTS OF FLOPPY DISK THAT ARE RELATED TO THE MSEDH 


A* \ 

>D1R 

' J. •.. . r :!n /e A r. a 3 . , < ]. ...\b *1 

MTF : M- F : F ■ < •••• ;-9- c.? » 1 

M3EDM OIF. - l - 09 -39 

7 i . « ■ ■' '? 79/. 90 by;t?'.v t- 


ms \M5EDM 
.> 0 1 R 



!. • 1 j . 1 

■ •' "> •• 

2 O’OdM 

,0 ! ,ibe l 






i ...7Q..UC 

1 i 

4 A . i 



L- 1 K .'■■■ 

1 --29-89 

1 t 

4 r J a 

LHl’UO 

c.; 

i i J. 

2-19-87 

10 

1 Oa 

MODSSUM 

FOR 

’-V !' 

2 - |9-S'7 

4 

29p 

MODSTRN 

0 9 

i £>47 

. 1 - ‘OR ■“ R 7 


49g 

MODSUM 

FOR 

4242 

2 -3 ~S / 


2-Va 

OAM 10 

E X E 

41394 

2 - ■ 1 ~ -87 


. 1 a 

MODS TRM 

F. ^ E 

42 733 

"t .■ \ i — ) -7 

\ 

'"Op 

i‘ ICO jLM 

E x E 

r ~ 3 3 A 2 

* T ... “7 


'■ -l 
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TABLE 3.3: EXAMPLE OF MSEDM 


4 


Computer Code 

Description 

Input Value 

MODSTRN 

Damage file used. 

DAM0UT3 . DD 


Problem title (= damage at N = 0.25 
and N = 0.5). 



Number of modes to use. 

12 


Number of elements for simply supported 
beam 

60 


Damage fraction. 

0.4 

OAMIO 

Damage file name. 

DAM0UT3 . DD 


Problem title (= damage at N * 0.25 
and N * 0.5) 



Number of modes to use. 

12 


Modal damage indicator values 
[scalars D(w,)], and corresponding 
critical values are listed below. 


MOOSUM 

Damage file used. 

DAM0UT3.DD 


Problem title (= damage at N * 0.25 
and N = 0.5) 



Number of modes to use. 

12 


Number of elements. 

60 


Damage fraction. 

0.4 


The output from the various computer codes are shown on the following pages: 
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TABLE 3.3 (continued) 


SIMPLY SUPPORTED BEAM 
DAMAGE FILE USED : DAM0UT3. DD 

DAMAGE TITLE s DAMAGE AT N = 0.25 AND AT N = 0.5 
NUMBER OF MODES t 12 
NUMBER OF ELEMENTS : 60 
DAMAGE FRACTION I .4000 


MODE 

MAXIMUM ELEMENT 

ELEMENT 

NO. 

STRAIN ENERGY 

NO. 

1 

1 . OOOE+OO 

31 

2 

1 . OOOE+OO 

45 

3 

1 . OOOE+OO 

31 

4 

1 . OOOE+OO 

38 

5 

1 . OOOE+OO 

31 

6 

1 . OOOE+OO 

45 

7 

1 . OOOE+OO 


8 

1 . OOOE+OO 

34 

9 

1 . OOOE+OO 

17 

10 

1 . OOOE+OO 

45 

1 1 

1 . OOOE+OO 

25 

12 

1 . OOOE+OO 

38 


DAMAGE 

FILENAME 

: DAM0UT3 . DD 

DAMAGE 

AT N = 0. 

25 AND 

AT N = 0 

NUMBER 

OF MODES 

: 12 


MODE 

D <W> 


Dcr l t ( W) 

1 

1 00 . 000 


10. 000 

dL 

100. OOO 


1 0 . 000 

“7 

1 00 . 000 


10. 000 

4 

1 . 000 


1 0. 000 

5 

1 00 . 000 


10 . 000 

6 

100. 000 


10. 000 

7 

100. 000 


10. 000 

8 

1 . 000 


1 0 . 000 

9 

IOC). 000 


1 0 . ooo 

10 

1 00 • 000 


1 0 . ooo 

1 1 

100. 000 


10. ooo 

12 

1 . 000 


1 0 . ooo 
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TABLE 3.3 (continued) 


SIMPLY SUPPORTED BEAM 
DAMAGE FILE USED t DAMOUT3. DD 

DAMAGE TITLE s DAMAGE AT N = 0.25 AND AT N = 0.5 
NUMBER OF MODES : 12 

NUMBER OF ELEMENTS : 60 
DAMAGE FRACTION : .4000 


ELEMENT DAMAGE SUMMARY 


ELEMENT 

NO. 


MODE 

NO. 


POSSIBLE 

DAMAGE 



1 

2 

T 

4 

er 

—1 

6 

7 

8 

9 

10 

1 1 

12 


1 

X 

X 

X 

A 

X 

X 

X 

X 

X 

X 

X 

X 

YES 


X 

X 

X 

X 

X 

y % 

X 

X 

X 

X 

X 


NO 

3 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 


NO 

4 

X 

X 

X 


X 

X 

X 


X 

X 

X 


NO 

3 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 

6 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 

7 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 


NO 

8 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 


NO 

9 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 


NO 

10 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 

1 1 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 

12 

X 

X 

X 


X 

X 

X 


X 

X 

X 


NO 

13 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 


NO 

14 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


NO 

13 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

YES 

16 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

YES 

17 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


NO 

18 

X 

X 

X 

X 

X 

X 

X 


X 

X 

X 


NO 

19 

X 

X 

X 


X 

X 

X 


X 

X 

X 


NO 

20 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 

21 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 

22 

:< 

X 

X 


X 

X 

X 

X 

X 

X 

X 


NO 

23 

X 

X 

X 


A 

K 

X 

X 

X 

X 

X 


NO 

24 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 


NO 

25 

X 

X 

X 


X 

X 

X 


X 

X 

X 

)■[ 

NO 

26 

x 

X 

X 


X 

X 

X 


X 

x 

X 

X 

NO 

27 

X 

X 

X 


X 

X 

X 


X 

X 

X 


NU 

28 

X 

X 

X 

X 

X 

X 

X 


A 




NO 

29 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 



MO 

30 

X 

A 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

YES 

31 

X 

A 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

YES 

32 

X 

A 

X 

X 

X 

X 

y 

X 

X 

X 

X 


Nu 

33 

X 

A 

X 

X 

X 

X 

X 


X 

X 

X 


NO 

34 

X 

X 

X 


X 

X 

;< 


X 

X 

M 


NO 

33 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 

36 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 

37 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 


NO 

38 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 


NO 

39 

X 

X 

X 


X 

X 

X 

X 

X 

X 

X 


NO 

40 

X 

X 

X 


X 

X 

X 


X 

X 

X 

X 

NO 


:s indicates possible modal element damage 
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TABLE 3.3 (continued) 


SIMPLY SUPPORTED BEAM 
DAMAGE FILE USED : DAM0UT3.DD 

DAMAGE TITLE : DAMAGE AT N = 0.25 AND AT N = 0.5 
NUMBER OF MODES : 12 

NUMBER OF ELEMENTS ! 60 

DAMAGE FRACTION s .4000 





ELEMENT 

DAMA 

GE 

SUMMARY 






ELEMENT 







MODE 





POSSIBLE 

NO, 







NO. 





DAriAGE 


1 

'**1 

"T 

4 

C T 

6 

7 a 

9 

10 

1 1 

12 


41 

X 

x 

X 


A 

x 

x 

X 

X 

A 

X 

r:o 

42 

X 

X 

X 


'A 

A 

X 

X 

A 

A 


NO 

43 

X 

X 

X 


A 

X 

X 

X 

X 

X 


NO 

44 

X 

X 

x 



X 

X X 

X 

X 

X 


NO 

45 

X 

X 

X 


A 

X 

X X 

X 

X 

X 

X 

YES 

46 

x 

X 

X 


A 

X 

X X 

X 

A 

X 

X 

YES 

47 

X 

X 

X 


5>! 

X 

X X 

X 

A 

X 


NO 

48 

X 

X 

X 

X 

X 

x 

X 

X 

A 

X 


NO 

49 

X 

X 

x 


)■( 

X 

X 

X 

X 

X 


NO 

50 

X 

X 

X 


X 

x 

X 

X 

X 

X 

X 

NO 

51 

X 

X 



X 

X 

X 

X 

X 

X 

X 

NO 

52 

X 

X 

x 


A 

X 

X X 

X 

X 

X 


NO 

53 

X 

X 

X 


X 

X 

X X 

A 

X 

X 


NO 

54 


X 

X 


X 

X 

X X 

X 

X 

X 


NO 

55 

X 

X 

X 


X 

X 

X 

X 

y t 

X 

X 

NO 

56 

X 

x 

X 


X 

X 

X 

X 

A 

X 

X 

NO 

57 

X 

X 

X 


X 

X 

X 

X 


X 


NO 

58 

X 

X 

X 

x 

X 

X 

X 

X 

A 

X 


NO 

59 

X 

X 

X 

X 

X 

X 

X X 

X 

A 

X 


NO 

60 

X 

X 

X 

X 

X 

X 

X X 

X 

X 

X 

X, 

YES 


x indicates possible modal element damage 
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4.0 CONTENTS OF FLOPPY DISKS 

A single floppy disk contains all the software necessary to solve 
either the module transfer function matrix*/damage assessment or modal 
strain energy distribution method problem defined in Reference 1. The 
software contained on the floppy is listed in Tables 2.2 and 3.2. Also, 

each floppy contains the input data for the two sample problems described 
earlier. The software on the floppies consists of the Fortran source 
listings and the corresponding executable codes. 

The floppies are 5-1/4-in., 1.2 MByte, high density 96 TPI. Three 

floppy disks, with software, have been provided with this report. 


* This is true provided that MAC/RAN IV has been installed on the micro- 
computer being used for these analyses. 
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APPENDIX A 

SOURCE CODE LISTINGS FOR 
MTFM AND MSEDM APPROACHES 





Contained in this appendix are the Fortran source code listings 
for the two damage assessment methods. 
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MODULE TRANSFER FUNCTION MATRIX (MTFM) 
SOURCE CODE LISTINGS 
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ooo oooooo ooo oo ooo 


PROGRAM NASA - PURPOSE : TO CREATE A BATCH FILE 
WHICH WILL EXECUTE MACRAN SEVERAL TIMES 


NASA. 


IMPLICIT INTEGER*2 (I-N) 

DIMENSION ICHNI (10), ICHN0( 14) 

CHARACTER*80 KEY 
CHARACTER* 30 TITLE 

CHARACTER* 12 MFIL1 I ,MFIL10,MFIL2I ,MFIL20 
CHARACTER* 12 MFIL3I,MFIL30,MFIL4I,MFIL40 
CHARACTER* 1 2 FILIN, FIL0T1 , FIL0T2 
CHARACTER*6 MFIL1 ,MFIL2 ,MFIL3,MFIL4,FIL0T 
CHARACTER*3 FX1 , FX2 , FX3 , FX4, AYN 
CHARACTER* 1 ATYPE, FILET( 12) 

L0GICAL*2 EXISTS, IER 

COMMON/NASADT/ICHNI , ICHNO, NINPUT , NOUTPT , FILIN , NPNTS 
EQUIVALENCE ( FILOT2 , FILET ( 1 ) ) 


CALL VCLEAR 
CALL VCURXY(O.O) 

DISPLAY PROGRAM NAME TO USER 

KEY = ’PROGRAM NASA - TO GENERATE A MACRAN BATCH FILE!’ 
CALL PRMPT( KEY, 1.00) 

BEGNIN HERE 


OPEN UP FILE CONTAINING NPTS 

OPEN(31,FILE=’ INTRFS.DAT’) 

READ(31 ,8831 ) NPTS 
8831 FORMAT (14) 

NPNTS = (NPTS+3)/4 
IR = 7 

10 CONTINUE 

KEY = ’ENTER INPUT MACRAN FILENAME (INCLUDING .SDF) : $’ 
CALL KEYBD( KEY, R, I,FILIN,3,(12),0,IR) 

SEE IF FILE EXISTS 

INQUIRE ( FI LE=FI LIN, EXIST: EXISTS) 

IF(EXISTS) GO TO 15 

IR = IR + 1 

CALL VCURXY(O.IR) 

WRITE(* ,9010) 

9010 FORMATOX, ’FILE DOES NOT EXISTS - PROGRAM ABORT!!’) 

STOP 

C 

C MAKE SURE IT HAS A .SDF EXTENSION 

C 

15 CONTINUE 

IR = IR + 1 
CALL VCURXY(0, IR) 

WRITE(*,8543) NPTS, NPNTS 
8543 FORMAT( IX, ’NUMBER OF POINTS : ’,14, 

1 /IX, ’NUMBER OF POINTS/SEGMENT ; ’,I4\) 

CALL VCURXY(P.IR) 

CALL FILECH( FILIN , IER) 

IF(IER) GO TO 20 
IR = IR + 2 
CALL VCURXY(O.IR) 

WRITE( * , 9002 ) 
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9002 FORMAT ( IX, ’NO . SDF EXTENSION - PROGRAM ABORT!!’) 

STOP 

C 

C HERE IF OK 
C 

20 CONTINUE 

IR = IR + 2 

KEY = ’ENTER OUTPUT FILENAME - NO FILE EXTENSIONS PLEASE :$ 
CALL KEYBO( KEY ,R, I , FI LOT ,3,(6),0,IR) 

C 

C VERIFY NO EXTENSION AND GET CURRENT LENGTH 
C 

IFLEN = 6 

CALL FILEX(FILOT, IFLEN) 

NFLEN = IFLEN + 6 
C 

C NOW GET INPUT/OUTPUT CHANNELS 
C 

30 CONTINUE 

IR = IR + 2 
CALL VCURXY(0, IR) 

KEY = ’ENTER THE NUMBER OF INPUT CHANNELS ( <= 10 ) :$ 

CALL KEYBD( KEY , R,NINPUT ,A,2,(2),0,IR) 

IF(NINPUT.LE. 10) GO TO 40 

IR = IR + 1 

CALL VCURXY(O.IR) 

WRITE( * , 9030) 

9030 FORMAT ( 1 X TOO MANY INPUT CHANNELS SPECIFIED’X) 

GO TO 30 
40 CONTINUE 

IR = IR + 2 
CALL VCURXY(O.IR) 

WRITE(* , 9031 ) 

9031 FORMATC IX, ’SPECIFY MACRAN CHANNEL TO USE FOR EACH INPUT \) 

IR = IR + 1 
CALL VCURXY(O.IR) 

KEY = ’ — > $’ 

DO 45 IC=1 .NINPUT 
IR = IR + 1 

CALL KEYBD(KEY,R, ICHNI (IC) ,A,2,(2),0,IR) 

45 CONTINUE 

C 

C NOW GET OUTPUT CHANNELS 

C 

48 CONTINUE 

CALL VCLEAR 
IR = 2 

CALL VCURXY(O.IR) 

KEY = ’ENTER THE NUMBER OF OUTPUT CHANNELS ( <= 14 ) :$ 
CALL KEYBD ( KEY , R , NOUTPT ,A,2,(2),0,IR) 

IF(NOUTPT . LE. 14) GO TO 50 

IR = IR + 1 

CALL VCURXY(O.IR) 

WRITE( * , 9040) 

9040 FORMAT ( 1 X ,’ TOO MANY OUTPUT CHANNELS SPECIFIED’X) 

GO TO 48 
C 

C NOW GET OUTPUT CHANNELS 

C 

50 CONTINUE 

IR = IR + 2 
CALL VCURXY(0 , IR) 

WRITE( * , 9041 ) 
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9041 FORMAT ( IX, 'SPECIFY MACRAN CHANNEL FOR EACH OUTPUT ’ \) 

IR = IR + 1 
CALL VCURXY(O.IR) 

1 KEY = $’ 

DO 55 IC = 1.N0UTPT 
IR = IR + 1 

CALL KEYBD( KEY,R, ICHN0( IC) ,A,2,(2),0,IR) 

55 CONTINUE 

C 

C SEE IF A DAMAGED OR UNDAMAGED CASE 

C 

IR = IR + 2 
CALL VCURXY(O.IR) 

54 CONTINUE 

KEY = 'IS THIS A (D)AMAGED OR (U)NDAMAGED RUN? :$’ 

CAL L KEYBD ( KEY , R , I , ATYPE ,3,(1),0,IR) 

IF(ATYPE. EQ. ’ D’ .OR. ATYPE. EQ. ’d’ ) GOTO 56 
IF(ATYPE. EQ. 'U' .OR. ATYPE. EQ. ’ u’ ) GO TO 56 
GO TO 54 

56 CONTINUE 

GET TITLE FOR THIS RUN 

KEY = 'ENTER A DESCRIPTIVE TITLE FOR THIS RUN : $’ 

IR = IR + 2 

CALL KEYBD (KEY , R , I , TITLE , 3 , ( 30) , 0 , IR) 

NOW START CREATING THE BATCH FILE 

OPEN ( 6 , FILE= ’ NASAD.BAT ’ , STATUS= ’ NEW’ ) 

OPEN ( 7 , FI LE= ’ NASAREN . BAT ’ , STATUS= ’ NEW ’ ) 

OPEN ( 8 , FI LE= ’ MATRX . DAT ’ , STATUS= ’ NEW ’ ) 

WRITE(8,4010) NINPUT.NOUTPT, ATYPE, TITLE 
4010 FORMAT ( 15 , 2X , I5,2X,A1/A30) 

ASSIGN GENERIC FILENAMES TO MACRAN FILES 

MFIL1 = 'TRNSD1 ’ 

MFIL2 = ’ TRNSD2 ’ 

MFIL3 = ’TRNSD3’ 

MFIL4 = 'TRNS04' 

FX1 = ’ INP’ 

FX2 = 'OUT' 

FX3 = 'SOF' 

FX4 = 'OAT' 

IL = 6 
C 

C NOW LOOP OVER THE NUMBER OF OUTPUT CHJANNELS TO PROCESS 
C 

DO 100 ICHN= 1 , NOUTPT 
C 

C GET FIRST MACRAN FILENAME 

C 

CALL FILEM(MFIL 1 ,MFIL1 I , FX1 , ICHN.IL) 

OPEN( 10 , FILE=MFIL1 I , STATUS= ’ NEW’ ) 

C 

C NOW WRITE OUT FIRST MACRAN CARD 

C 

CALL MACRD1 (FILIN , NINPUT , ICHNI , ICHNO(ICHN) .TITLE , NPNTS) 
CLOSE ( 10 , STATUS: ’ KEEP’ ) 

CALL FILEM(MFIL1 , MFIL10.FX2 , ICHN, IL) 

C 
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C GET SECOND MACRAN FILENAME 

c 

CALL FILEM(MFIL2,MFIL2I, FXI.ICHN.IL) 
i OPEN( 10 , FILE=MFIL2I , STATUS= ’ NEW’ ) 

C 

C WRITE OUT SECOND CARD 
C 

CALL MACRD2 

CLOSE ( 1 0 , STATUS= ’ KEEP ’ ) 

CALL FILEM(MFIL2,MFIL20, FX2 , ICHN, IL) 

C 

C GET THIRD MACRAN FILENAME 

C 

CALL FILEM(MFIL3 , MFIL3I , FX1 , ICHN, IL) 
0PEN( 10 , FILE=MFIL3I ,STATUS= ’ NEW’ ) 

WRITE OUT THIRD MACRAN CARD 

CALL MACRD3 

CLOSE ( 1 0 , STATUS= ’ KEEP’ ) 

CALL FILEM(MFIL3,MFIL30,FX2, ICHN , IL) 

NOW ADD TO THE NASAD.BAT FILE 

WRITE(6 , 6010) MFIL1I.MFIL10.MFIL10 
WRITE(6 , 6060) MFIL2I .MFIL20.MFIL20 
WRITE( 6 ,6010) MFIL3I .MFIL30 , MFIL30 

GET OUTPUT FILENAME 

CALL FI LEM ( FI LOT , FI LOT 1 , FX3 , ICHN , IFLEN) 
CALL FILEM(FIL0T,FIL0T2,FX4, ICHN, IFLEN) 

WRITE OUT BATCH COMMANDS 

WRITE(6 , 6020) FI LOT 1 , FI LOT 1 
WRITE(6 , 6030) FILOT2 , FILOT2 
WRITE(8 , 4020 ) ( FILET(IFM) , IFM= 1 , NFLEN) 

100 CONTINUE 

4020 FORMAT (12A1) 

IF(IFLEN.EQ.6) GO TO 150 
CALL FILADCFILOT, IFLEN) 

150 CONTINUE 

FINISH OFF BATCH COMMANDS FOR NASA. BAT 

WRITE(6 , 6070) 

CLOSE(6,STATUS= ’KEEP’ ) 

NOW ADD TO NASAREN.BAT 

WRITE(7 , 6050) FILOT , FILOT 
CLOSE(7, STATUS= ’KEEP’) 

C 

C END IT ALL 
C 

CLOSE ( 7 , STATUS: ’ KEEP’ ) 

CALL VCURXY(0,23) 

STOP 
C 
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FORMAT SECTION 

6010 FORMAT (’MAC423 /R 10000 < * , A 1 2 , ’ >0: ’ ,A12/’DEl D:\A12) 
6020 FORMAT ( ’DEL’ , IX, A 12/ ’REN TRNSMD2.SDF’ , 1X,A12) 

6030 FORMAT ( ’DEL’ , IX, A 12/ ’REN HD2 . DAT ’ , IX, A12) 

6040 FORMAT ( ’ DEL ’ , 1 X , A6 , ’ 0 1 . BFT ’ 

1 /’REN FORDI.DAT’, IX, A6, ’01. BFT’) 

6045 FORMAT ( ’ DEL ’ , 1 X , A6 , ’ 0 1 . DFT ’ 

1 /’REN FORD0 1 . DAT ’ , 1 X , A6 , ’ 0 1 . DFT ’ ) 

6046 FORMAT ( ’ DEL ’ , IX, A6 , ’ 02 . DFT ’ 

1 /’REN FORD02.DAT’, IX, A6, ’02. DFT’ ) 

6050 FORMAT ( ’ DEL ’ , 1 X , A6 , ’ . BAT ’ 

1 /’REN NASAD.BAT’, IX, A6, ’ .BAT’ ) 

6060 FORMAT ( ’MACTRN /R 10000 <’,A12,’ > D:’,A12/’DEL D:’,A12) 
6070 FORMAT ( ’ TESTRD . EXE ’ / ’ NASAREN . BAT ’ ) 

END 
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SUBROUTINE MACCRD3 - PURPOSE TO WRITE OUT 3RD MACRAN CARD 


SUBROUTINE MACRD3 
IMPLICIT INTEGER*2 (I-N) 
WRITE( 10,9010) 

9010 FORMAT (’ FILIN TRNSMD2 . SDF’ 

1 /’FI LOUT HD2 . DAT ’ 

2 / ’ CTAPE ’ , 52X, 1 ’ 

3 /’END’ 

4 /) 

RETURN 

END 


Mftc£D3.Fe£_ 
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SUBROUTINE MACCRD2 - PURPOSE TO WRITE OUT 2ND MACRAN CARD 


SUBROUTINE MACRD2 
IMPLICIT INTEGER*2 (I-N) 

WRITE OUT STUFF 


Macros. Fc*l 


WRITE( 10,9010) 

9010 FORMAT( ’FILIN MATRXD2.SDF’ 

1 /’FI LOUT TRNSMD2.SDF’ 

2 /’TRANS TAPE’ 

3 /’END’ 

4 /) 

RETURN 

END 
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C 

C SUBROUTINE MACCR01 - PURPOSE TO WRITE OUT FIRST MACRAN CARD 
C 

SUBROUTINE MACRO 1 (FILIN , NINPUT , ICHNI , ICHNO, TITLE , NPNTS) 
IMPLICIT INTEGER*2 (I-N) 

DIMENSION ICHNI ( 10) 

CHARACTER* 12 FILIN 
CHARACTER*30 TITLE 
C 

C WRITE OUT FIRST COUPLE OF COMMENT CARDS 
C 

WRITEC 10,9010) TITLE 

9010 FORMAT (’JBHEAD ’ ,A30, 5X, ’TRANSFER FUNCTIONS’) 

WRITE( 10 , 9020) FILIN 
9020 FORMAT( ’FILIN’ ,2X,A12) 

C 

C NOW WRITE OUT CHANNEL NUMBERS 
C 

IFCNINPUT.GT.2) GO TO 50 

WRITEC 10, 9030) ICHNO, ( ICHNI ( IL) , IL=1 , NINPUT) 

GO TO 100 
50 CONTINUE 

WRITEC 10, 9030) ICHNO, (ICHNI (IL) , IL=1 ,2) 

NLINES = (NINPUT+1 )/4 
DO 60 IJ=1, NLINES 
IL = 3 + 4 * (IJ-1) 

IK = IL + 3 

IFCIK.GT. NINPUT) IK = NINPUT 
WRITE(1 0,9040) ( ICHNI ( IM) , IM=IL , IK) 

60 CONTINUE 

9030 FORMAT (’FILIN’ , IX , ’CCHAN’ ,14X,’1. ’ , IX, 3(9X, 12 , ’ . ’)) 

9040 FORMAT ( ’ FI L IN ’ , 1 X , ’ CCHAN ’ , 5X , 4( 9X , I 2 , ’ . ’ ) ) 

C 

C HERE FOR REST OF DATA 

C 

100 CONTINUE 

WRITE( 1 0 , 9050) NPNTS 
9050 FORMAT ( 

1 ’FILOUT MATRXD2 . SDF’ 

2/ ’ PSD TRANS’, 1 IX, 14,’.’, 

3/ ’PSD TOUTMTFC’ 

4/ ’PSD RETAIN’ 

5/ ’END’ 

6 /) 

RETURN 

END 
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SUBROUTINE FILMAP - TO CONSTRUCT A FILENAME 
CALLED FORM FILEM 

SUBROUTINE FILMAP(FIL1 , FIL2 , FX.ANUM, ILEN) 
IMPLICIT INTEGER*2 (I-N) 

CHARACTER*! ANUM(2) , FX(3) , FIL1 ( 12) ,FIL2( 12) 


DO 10 I=1,ILEN 
FIL2CI) = FILI(I) 

10 CONTINUE 

DO 15 1=1,2 

FIL2CILEN+I) = ANUM(I) 

15 CONTINUE 

FIL2CILEN+3) = ’ . ’ 

DO 20 1=1,3 

FIL2( ILEN+3+I ) = FX(I) 

20 CONTINUE 

CHANGE BLANKS TO ZEROES 

NLEN = ILEN + 6 
DO 30 1 = 1, NLEN 

IF(FIL2(I) . EQ. ’ ’) FIL2(I) = ’O’ 
30 CONTINUE 

RETURN 
END 
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SUBROUTINE TO CHECK ON FILE AND ITS LENGTH P| Lt/. FC 

SUBROUTINE FI LEX( FIL 1 , ILEN) 

IMPLICIT INTEGER*2 (I-N) 

CHARACTER* 1 FIL1(20) 

GET LENGTH OF FILE 

DO 5 1=1, ILEN 

I F ( FI L 1 (I) . EQ. ’ ’ ) GO TO 7 
5 CONTINUE 

GO TO 8 

7 CONTINUE 
ILEN = 1-1 

NOW SEE IF AN EXTENSION WAS ADDED 

8 CONTINUE 

DO 10 1=1, ILEN 

I F ( FI L 1 (I) . EQ. ’ . ’ ) GO TO 20 
10 CONTINUE 

RETURN 

20 CONTINUE 

ILEN = ILEN-1 

RETURN 

END 
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SUBROUTINE FILEM - TO CONSTRUCT A FILENAME F ( L£E“M . 

SUBROUTINE FILEM(FIL1 , FIL2 , FX, INUM, ILEN) 

IMPLICIT INTEGER*2 (I-N) 

CHARACTER* 1 FIL 1 ( 1 2 ) , FIL2 ( 1 2) , FX(3) 

CHARACTER*2 ANUM 

WRITE OUT INUM TO ANUM 

WRITE(ANUM, 900) INUM 
900 FORMAT (12) 

CALL FILCHG TO GET FILENAME 

CALL FILMAP(FIL1 ,FIL2,FX, ANUM, ILEN) 

RETURN 

ENO 
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SUBROUTINE FILECH - PURPOSE TO SEE IF A . SDF EXTENSION IS ON 

SUBROUTINE FI L ECH ( FI L 1 , IER) 

IMPLICIT INTEGERS (I-N) 

CHARACTER* 1 FIL 1 ( 1 2 ) 

LOGICAL*2 IER 

FIRST LOOK FOR A . 

IER = .FALSE. 

DO 10 1=1,12 

IFCFILI(I) . EQ. ’ . ’ ) GO TO 20 
10 CONTINUE 

RETURN 

20 CONTINUE 

ILEN = I 

I F ( FI L 1 (ILEN+1 ) .NE. ’S’ .AND. FIL1 (ILEN+1 ) . NE. ’ s’ )RETURN 
I F ( FI L 1 (ILEN+2) .NE. ’ D’ .AND. FI LI (ILEN+2) .NE. ’d’ ) RETURN 
I F ( FI L 1 (ILEN+3) .NE. ’F’ .AND. FIL 1 ( ILEN+3) . NE. ' f ’ )RETURN 
IER = .TRUE. 

RETURN 

END 


F/lEch-. kor. 
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INOFLOATCALIS 

CHARACTER*80 KEY 
CHARACTER* 14 FLNAME 


CALL VCLEAR 
CALL VCURXY(0 , 0 ) 

KEY = ’ENTER FILENAME TO DISPLAY : $’ 
CALL KEYBDC KEY , R , I , FLNAME ,3, (14), 0,5) 
OPEN( 7 , FILE=FLNAME , FORM= ’ UNFORMATTED’ ) 
10 CONTINUE 

READ(7) FTR, FTI 
WRITE(*,8000) FTR, FTI 
8000 FORMAT ( 2X , 1 PE 1 0 . 3 , 5X , 1 PE 1 0 . 3 ) 

GO TO 10 

STOP 

END 


Fl 
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SUBROUTINE FILAD( FILOT , IFLEN) F| L. D. ^Of<L 

SUBROUTINE FILAD(FILOT , IFLEN) 

IMPLICIT INTEGER*2 (I-N) 

CHARACTER* 1 FILOT(6) 

IFL = IFLEN + 1 
DO 10 I = I FL , 6 
FILOT ( I ) = ’0’ 

10 CONTINUE 

RETURN 
END 
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JNOFLOATCALLS . . - . c 

SUBROUTINE MACWRT (DAM , NFREQ, DELF .TITLE) MftC WK.T. r~CK. 

IMPLICIT INTEGERS ( I-N) 

COMPLEX DAM( 7 , 2500) , XDAM( 7 ) 

REALM RDAM( 16) 

CHARACTERM ADAM(7) 

CHARACTER*30 TITLE 

SET UP CHARACTER STRINGS 

ADAM(1) = ’ 0AM 1 ’ 

ADAM ( 2 ) = ’ DAM2 ’ 

ADAM(3) = ’ DAM3’ 

ADAM(4) = ’ DAM4 ’ 

ADAM( 5 ) = ’ DAM5 ’ 

ADAM(6) = ’ DAM6 ’ 

ADAM ( 7 ) = ’ DAM7 ’ 

OPEN UP OUTPUT DAMAGE. SCF FILE 

OPEN( 20 ,FILE=’ DAMAGE . SCF’ , STATUS= ’ NEW’ ) 

START WRITING OUT TO FILE 

WRITE(20,9090) TITLE 
9090 FORMAT (’JBHEAD’, IX, ’COMPUTED DAMAGE - \A30) 

WRITEC20 , 9050) 

WRITE( 20 , 9000) 

WRITE(20,9010) 

WRITE(20 , 9020) 

DO 20 1=1,6 
DO 10 IW= 1 , 2 

WRITE(20,9030) ADAM ( I ) 

10 CONTINUE 

20 CONTINUE 

DO 30 1=1,16 
RDAM(I) =0.0 
30 CONTINUE 

CALL GCURXY(IC.IR) 

IR = IR + 1 
CALL VCURXY(O.IR) 

DO 100 IFQ=1 , NFREQ 

RDAM(1) = FLOAT ( IFQ-1 ) * DELF 
DO 40 1=1,6 

XDAM(I) = DAM ( I , I FQ ) 

RI = AIMAG(XDAM(I)) 

RR = REAL(XDAM(I) ) 

RDAM(I*2 ) = RR 
RDAM(I*2+1) = RI 
40 CONTINUE 

WRITE(20,9040) (RDAM(I) ,1=1,16) 

100 CONTINUE 

WRITE(20 , 9060) 

CLOSE (20 , STATUS= ’KEEP’ ) 

9000 FORMAT ( ’CTAPE’ ,2X, ’ NASA 1311.05’) 

9010 FORMAT( ’ CTAPE ’ ,2X, 'DAMAGE ’ ,16X,’13’) 

9020 FORMAT ( ’ CTAPE ’ , 2X , ’ FREQUENCY ’ ) 

9030 FORMAT ( ’CTAPE’ ,2X,A4) 

9040 FORMAT ( ’CTAPE’ , 19X,4(E11 .3, IX) , 

1 /’CTAPE C’ ,17X,4(E11.3,1X), 

2 /’CTAPE C’ , 1 7X , 4(E1 1 . 3 , IX) , 

3 /’CTAPE C ’ , 1 7X , 4(E1 1 . 3 , IX) ) 

9050 FORMATC ’FI LOUT DAMAGE. SDF’ ) 

9060 FORMAT (’ END ’ ) 

RETURN 
END 
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SUBROUTINE DAMCOM(IFQ) DA-MCGM Ffl 

IMPLICIT INTEGER*2 (I-N) 

COMPLEX*8 TFU (14,10), TFD(14,10) 

COMPL EX*8 0AM ( 7 , 2 500 ) , AUX 1 , AUX2 , AUX3 , AUX4 , AUX5 
DIMENSION ILMI (14) , ILMO( 14) 

COMMON/DAMC/T FU , T FD , DAM , I LM I , I LMO , N I NPUT , NOUT PT 

ZERO OUT AUX 

AUX1 = (0. ,0. ) 

AUX2 = (0.,0.) 

AUX3 = (0.,0.) 

AUX4 = (0. ,0. ) 

AUX5 = (0. ,0. ) 

COMPUTE DAMAGE 

DO 20 KK= 1 ,NOUTPT 
DO 19 IK= 1 , NOUTPT 

IF(KK. EQ. ILMO( IK) ) GO TO 20 

19 CONTINUE 

DO 10 JJ= 1 , NINPUT 
DO 11 IJ=1, NINPUT 

IF(JJ. EQ. ILMI(IJ) ) GO TO 10 
11 CONTINUE 

C AUX1 = AUX1 + TFU(KK, JJ) * FTO(KK) * FTI(JJ) 

C AUX2 = AUX2 + TFD(KK.JJ) * FTO(KK) * FTI(JJ) 

AUX3 = AUX3 + TFU(KK, JJ) * TFU(KK.JJ) 

AUX4 = AUX4 + TFU(KK, JJ) * TFD(KK.JJ) 

AUX5 = AUX5 + TFD(KK, JJ) * TFD(KK.JJ) 

10 CONTINUE 

20 CONTINUE 

DAM ( 1 , I FQ ) = (AUX5 + AUX4) / (AUX4 + AUX3) 

DAM(2, IFQ) = (AUX4 + AUX3) / (AUX5 + AUX4) 

DAM ( 3 , I FQ ) = AUX3/AUX4 
DAM(4, IFQ) = AUX4 / AUX 3 
DAM(5 , IFQ) = AUX4 / AUX5 
DAM(6, IFQ) = AUX5 / AUX4 
RETURN 
END 
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JNOFLOATCALLS 

IMPLICIT INTEGER*2 (I-N) 
C0MPLEX*8 TFU(14,10),TFD (14,10) 
C0MPLEX*8 DAM(7 , 2500) 

DIMENSION I LM 1(14) , ILMO( 14) 
CHARACTER*80 KEY 
CHARACTER*30 TITLE 
CHARACTER*20 TFUF.TFDF 
CHARACTER*! ATYPE 


Damage, fqi^ 


COMMON/DAMC/TFU , TFD , DAM , I LMI , I LMO, NINPUT , NOUTPT 


CALL VCLEAR 
CALL VCURXY(0,0) 

IR = 0 

GET EACH FILENAME, OPEN IT UP AND READ HEADER INFO 

TFUF = ’HWUMAT.DAT’ 

OPEN(6,FILE=TFUF,FORM=’ UNFORMATTED’ ) 

READ ( 6 ) TITLE.NFREQ, NINPUT, NOUTPT , ATYPE , DELF 
WRITE(*,9000) TITLE, NINPUT, NOUTPT, NFREQ, DELF 
9000 FORMAT ( IX, A30 

1 /IX, ’NUMBER OF INPUT CHANNELS 

2 /IX, 'NUMBER OF OUTPUT CHANNELS 

3 /IX, ’NUMBER OF FREQUENCY POINTS 

4 /IX, ’FREQUENCY STEP SIZE 

IR = IR + 6 
CALL VCURXY(0, IR) 

IF(ATYPE. EQ. ’U’ .OR. ATYPE. EQ. ’u’ ) 

IF(ATYPE.EQ. ’D’ .OR. ATYPE. EQ. ’ d ’ ) 

9010 FORMAT (IX,’ UNDAMAGED ’ \ ) 

9015 FORMAT ( IX, ’DAMAGED’ \) 

IR = IR + 2 
CALL VCURXY(O.IR) 

C 

C GET DAMAGED TF FILENAME 

C 

TFOF = ’HWDMAT.DAT’ 

OPEN (7 , FILE=TFDF, FORM= ’ UNFORMATTED’ ) 

READ(7 ) TITLE , NFREQ, NINPUT , NOUTPT .ATYPE , DELF 
WRITE(* , 9000) TITLE, NINPUT, NOUTPT, NFREQ, DELF 
IR = IR + 6 
CALL VCURXY(O.IR) 

IF(ATYPE. EQ. ’U* .OR. ATYPE. EQ. ’ u ’ ) WRITE(*,9010) 
IF(ATYPE. EQ. ’ D’ .OR. ATYPE. EQ. ’d’ ) WRITE(*,9015) 
IR = IR + 2 

C 

C SEE ABOUT ROW/COLUMN ELIMINATION 

C 

PAUSE 

CALL VCLEAR 
IR = 0 

CALL VCURXY(O.IR) 

DO 40 1=1,14 
ILMI(I) = 0 
ILMO(I) = 0 
40 CONTINUE 


\I4 

\I4 

\I4 

’ , 1PE10.3, ’ HZ’\) 


WRITE(* , 9010) 
WRITE(*, 9015) 
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KEY = 'HOW MANY INPUT CHANNELS DO YOU WANT TO ELIMINATE? 
CALL KEYBD( KEY ,R, IELMI ,A,2,(2),0,IR) 

IR = IR + 1 

IF(IELMI.EQ.O) GO TO 75 
CALL VCURXY(O.IR) 

WRITE(* , 9200) 

9200 FORMAT (IX,’ ENTER CHANNELS TO ELIMI^ATE’X) 

IR = IR + 1 
CALL VCURXY(O.IR) 

DO 50 1=1, IELMI 
KEY = $’ 

CALL KEYB0(KEY ,R , ILMI (I),A,2,(2),0,IR) 

IR = IR + 1 
50 CONTINUE 

75 CONTINUE 

IR = IR + 1 

KEY= ' HOW MANY OUTPUT CAHNNELS DO YOU WANT TO ELIMINATE? 
CALL KEYBD(KEY , R , I ELMO , A, 2 , (2) , 0 , IR) 

IF(IELMO.EQ.O) GO TO 100 

IR = IR + 1 

CALL VCURXY(0, IR) 

WRITE(* , 9200) 

IR = IR + 1 
CALL VCURXY(O.IR) 

DO 80 1=1 .IELMO 
KEY = ’~> $’ 

CALL KEYBD(KEY,R, ILMO(I) , A, 2 , (2) ,0, IR) 

IR = IR + 1 
80 CONTINUE 

100 CONTINUE 

NOW READ IN ALL THE DATA, ONE FREQ POINT AT A TIME AND COMPUTE 
DAMAGE FOR THAT FREQ POINT 

NTOT = NINPUT + NOUTPT 
DO 150 I FQ= 1 ,NFREQ 
DO 110 10=1, NOUTPT 
READC6) (TFU(IO, I), 1=1, NINPUT) 

REA0(7) (TFDCIO, I), 1=1, NINPUT) 

110 CONTINUE 

CALL DAMCOM(IFQ) 

150 CONTINUE 

CALL MACWRT TO WRITE OUT DATA 

CALL MACWRT ( DAM , N FREQ , DE L F , TITLE) 

END IT ALL 

CALL GCURXY(IC.IR) 

IR = IR + 1 
CALL VCURXY(0 , IR) 

STOP 
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INOFLOATCALLS 

IMPLICIT INTEGER*2 (I-N) 

COMPLEX*8 HW(30) , FT(30) 

DIMENSION ADAT ( 30 ) , ODAT 1(30), ODAT 2(30), FREQD( 2 ) 
CHARACTER* 12 FLNAME 
CHARACTER* 1 ADUM.ATYPE 
CHARACTER*30 TITLE 




C READ IN DATA FROM FILE 

C 

0PEN(8,FILE=’MATRX.DAT’ ) 

READ(8 , 5000) NINPUT ,NOUTPT , ATYPE .TITLE 
5000 FORMAT ( 15 , 2X, 15 , 2X , A1/A30) 


C 

C NOW READ IN EACH FILENAME AND OPEN UP A CHANNEL 
C 

DO 2 1=1 ,NOUTPT 

READ(8, 5010) FLNAME 
ICHN =1+20 
OPEN(ICHN, FILE=FLNAME) 

2 CONTINUE 

5010 FORMAT (A1 2) 

CLOSE (8, STATUS= ’KEEP’) 


C 

C OPEN UP 
C 


4 

5 


OUTPUT FILE 


IF(ATYPE. EQ. ’ U’ .OR. ATYPE. EQ. ’ u’ ) GOTO 
OPEN (6 , FILE= ’ HWDMAT . DAT ’ , STATUS: ’ NEW’ , FORM= ’UNFORMATTED 


GO TO 5 

OPEN (6, FILE= ’HWUMAT.DAT’ , STATUS: ’NEW’ , FORM= ’ UNFORMATTED’ ) 


CONTINUE 


C READ IN FIRST RECORD AND GET THE NUMBER OF FREQ RECORDS 
C 

READ(21 ,9000) NREC 
9000 FORMAT (/30X, 15) 

NFREC = (NREC-1 )/2 
NRECD = NREC + 2 
REWINO 21 

NLINE = (NREC+3)/4 


C 

C READ CHANNEL 21 AND COMPUTE 4 OF FREQ POINTS 


DO 11 1=1, NRECD 
READ(21 ,9001) ADUM 


11 

CONTINUE 
NFREQ = 0 


12 

CONTINUE 

DO 13 1=1, NLINE 



READ(21 ,9001 ,END=14) 

ADUM 

13 

CONTINUE 

NFREQ = NFREQ + 1 
GO TO 12 


14 

CONTINUE 
REWIND 21 

WRITE(*,4010) NFREQ 

POINTS READ’) 

4010 

FORMAT (IX, I 5, ’ FREQUENCY 
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C READ IN FIRST TWO FREQ VALUES AND COMPUTED DELTA FREQ 
C 

DO 21 1=1 ,NRECD 
READ(21 ,9001) ADUM 

21 CONTINUE 

DO 22 1=1,2 

READ(21 ,5800) FREQD(I) 

DO 22 I L = 2 , NLINE 
READC21 ,9001) ADUM 

22 CONTINUE 

5800 FORMAT (24X , El 1 . 3) 

REWINO 21 

DELF = FREQD(2) - FREQD(I) 

WRITE(*,5810) DELF 

5810 FORMAT( IX, ’FREQUENCY STEP SIZE : \1PE10.3,’ HZ’) 

C 

C READ OVER DUMMY DATA IN EACH FILE 
C 

DO 6 IFL= 1 , NOUTPT 
ICHN = 20 + IFL 
DO 10 1=1 , NRECD 

READ( ICHN, 9001) ADUM 
10 CONTINUE 

6 CONTINUE 

9001 FORMAT (A1) 

C 

C NOW WRITE OUT HEADER INFO 
C 

WRITE(6) TITLE, NFREQ.NINPUT, NOUTPT, ATYPE, DELF 
C 

C NOW LOOP OVER EACH FILE, READING A ROW AT EACH FREQ POINT 
C 

20 CONTINUE 

DO 30 I FL=1, NOUTPT 
ICHN = IFL + 20 
DO 15 IW=1, NLINE 
II = ( I W— 1 ) *4+1 
JJ = II + 3 

READ( ICHN, 901 0,END= 100) (ADAT(K) , K=II , JJ) 

15 CONTINUE 

C 

C NOW PLACE IN COMPLEX ARRARY 
C 

DO 18 IW= 1 , NFREC 

HW(IW) = CMPLX(ADAT(IW*2) , ADAT ( IW*2+ 1 ) ) 

18 CONTINUE 

WRITE(6) (HW(I), 1=1, NFREC) 

30 CONTINUE 

9010 FORMAT (24X,4(E1 1 . 3, IX)) 

GO TO 20 
100 CONTINUE 

DO 105 IFL=1, NOUTPT 
ICHN = 20 + IFL 
CLOSE ( ICHN , STATUS= ’ KEEP’ ) 

105 CONTINUE 

STOP 
END 
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INOFLOATCALLS 

SUBROUTINE OUT1 OUT J.FO/C 

i PARAMETER ( NT IME= 5000 ) 

IMPLICIT INTEGER*2 (I-N) 

DIMENSION XIN(24,NTIME) 

CHARACTER* 12 MACOUT , MACBI N , MAC IN 
CHARACTER*8 CHANDES(24) 

COMMON/NASA/XIN , CHANDES , NTOT , NPTS , DT ,MACIN, MACOUT , MACBIN 
NREC = NTOT + 1 
C 

C WRITE OUT FIRST SET OF CTAPE CARDS 
C 

0PEN(4, FILE=MACOUT , STATUS= ’NEW’ ) 

WRITE(4, 9001 ) MACIN 

9001 FORMAT (’JBHEAD’ , IX, ’TIME HISTORIES FROM COSMOS FILE - \A12) 

WRITE(4, 9000) MACBIN 
9000 FORMAT (’FI LOUT ’ , A12 ) 

WRITE(4, 9010) 

9010 FORMAT ( ’CTAPE’ ,2X, ’NASA 1311.05’) 

WRITE(4, 9020) NREC 

9020 FORMAT ( ’ CTAPE ’ , 2X , ’ COSMOS/M ’ , 1 5X , 12 , ’ . ’ ) 

WRITE(4, 9030) 

9030 FORMAT ( ’CTAPE’ ,2X, ’TIME’) 

DO 10 1=1, NTOT 

WRITE(4,9040) CHANDES ( I ) 

10 CONTINUE 

9040 FORMAT (’ CTAPE \2X, AS) 

C 

C WRITE OUT DATA 
C 

IF(NREC. GT. 4) GO TO 200 
C 

C HERE IF LESS THAN 4 OUTPUT CHANNELS 
C 

DO 110 1=1, NPTS 
TIME = (1-1) * DT 

WRITE(4,9050) TIME, (XIN( IFN, I ) , IFN= 1 ,NTOT) 

110 CONTINUE 

WRITE(4, 9090) 

RETURN 

200 CONTINUE 

NREC = NREC-4 
NLINE = (NREC+3)/4 
DO 210 1=1, NPTS 
TIME = (1-1) * DT 

WRITE(4, 9050) TIME, (XIN( IFN , I ) , IFN= 1 , 3) 

DO 215 IP = 1, NLINE 
II = (IP-1) *4+4 
JJ = II + 3 

IF(JJ.GT.NTOT) JJ = NTOT 

WRITE(4, 9060) (XIN( IFN , I ) , IFN=II , JJ ) 

215 CONTINUE 

210 CONTINUE 

WRITE(4,9090) 

9090 FORMAT ( ’ END’ ) 

9050 FORMAT( ’ CTAPE’, 19X, 4(1 PE1 1.3, IX)) 

9060 FORMAT ( ’CTAPE C ’ , 1 7X , 4( 1 PE1 1 . 3 , IX) ) 

RETURN 

END 
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C THIS PROGRAM TAKES THE OUTPUT DATA (*.PLT) FILE FROM COSMOS * iNT&FS.K^ 

C AND PREPARES AN INPUT FILE (*.SCF) FOR MACRAN. * 

i C * 

C NASA PROJECT JOB » 1311.05B * 

C * 

c * 

c * 

C THE PROGRAM WILL OPEN A MAXIMUM OF 25 CHANNELS. TIME WILL * 

C ALWAYS BE CHANNEL 1 FOLLOWED BY 10 INPUT CHANNELS AND THE * 

C OUTPUT CHANNELS ( MAXIMUM = 14 ). THE PROGRAM ALSO CREATES A * 

C BATCH FILE NAMED NASAMAC.BAT, WHICH UPON EXECUTION, BESIDES * 

C PERFORMING THE AFOREMENTIONED TASK WILL CALL MACRAN AND * 

C CONVERT THE ASCII(*.SCF) FILE TO A BINARY(*.SDF) FILE. * 


* 

* 

* 

DEFINITION OF VARIABLES * 

* 

TIME : DISCRETE TIME SERIES * 

XLOD : INPUTS FROM COSMOS * 

XOUT : OUTPUTS FROM COSMOS * 


C NTIME : PARAMETER FOR VARIABLE DIMENSION, * 

C : SET TO MAXIMUM OF 5000 POINTS. * 

C NPTS : NO OF POINTS IN THE TIME HISTORY * 

C ICHANO : NO. OF OUTPUT CHANNELS * 

C ICHANI : NO. OF INPUT CHANNELS * 

C * 

C * 

C * 

C USER INPUTS * 

C * 

C NAME OF THE INPUT FILE NAME : (*.PLT) * 

C NAME OF THE MACRAN ASH FILE : (*.SCF) * 

C NAME OF THE MACRAN BINARY FILE : (*.SDF) * 

C NO OF POINTS IN THE TIME HISTORY : NPTS * 

C * 

C * 

C * 

C NOTES * 

C * 

C 1. THE SOURCE CODE HAS BEEN PROGRAMMED TO ACCEPT 5000 TIME * 

C POINTS, MADE AVAILABLE TO IT FROM COSMOS PLOT FILE * 

C O.PLT). SHOULD IT BE NECESSARY TO INCREASE THIS LIMIT * 

C CHANGE THE VALUE OF ’NTIME= ’ APPEARING IN THE PARAMETER * 

C STATETMENT OF THE MAIN PROGRAM — REFORM, AND THE FIVE * 

C SUBROUTINES IT CALLS NAMELY; OUT1, OUT2, OUT3, 0UT4, * 

C AND OUT 5. THIS HAS TO BE DONE BY EDITING THE SOURCE * 

C PROGRAM INTRFS. FOR . * 

C * 

C 2. ITS IMPERATIVE THAT THE USER ACQUAINTS HIMSELF WITH THE * 

C COSMOS PLOT FILE AND IS FAMILIAR WITH THE INPUTS AND THE * 

C OUTPUTS SPECIFIED IN THE COSMOS RUN. THE CODE IS WRITTEN * 

C FOR TEN INPUTS(STRESSES, DISPLACEMENTS). IT CONSIDERS * 

C ACCELERATION AS ITS OUTPUTS ONLY. IT IS RECOMMENDED * 

C THAT THE USER MAKES A NOTE OF THESE PARAMETERS AND * 

C UTILIZES ONLY THE NON-ZERO QUANTITIES IN THE SECOND MACRO. * 

C THIS WILL PRECLUDE ERRONEOUS TRANSFER FUNCTIONS RESULTING * 

C FROM DIVISION OR MULTIPLICATION OF ZERO. * 

C * 

C 3. TO EXECUTE THIS PROGRAM RUN NASAMAC.BAT * 

C * 


C********** *************************************************** ******* 
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$NOFLOATCALLS 

PARAMETER(NTIME=5000) 

IMPLICIT INTEGER*2 (I-N) 

CHARACTER* 120 ALINE 
CHARACTER* 12 MAC IN , MACOUT , MACBIN 
CHARACTER*4 ATIME 
CHARACTER* 3 AYN 
CHARACTER*80 KEY 
DIMENSION XIN(24 , NTIME) 

CHARACTER*8 CHANDES(24) 

COMMON/NASA/XIN , CHANDES , NTOT , NPTS , DT ,MACIN, MACOUT , MACBIN 
C 

C CLEAR SCREEN, OISPLAY PROMPT ,ETC 
C 

CALL VCLEAR 
CALL VCURXY(0,0) 

KEY= ’ INTRFS - PROGRAM TO REFORM COSMOS FILE TO MACRAN FILE $ 
CALL PRMPTCKEY, 1 .00) 

CALL GCURXY(IC.IR) 

IR = IR + 1 
C 

C GET INPUT AND OUTPUT FILENAMES 

Q 

KEY = ’ENTER COSMOS PLOT FILENAME TO PROCESS : $’ 

CALL KEYBDCKEY , R , I ,MACIN ,3,(12),0,IR) 

IR = IR + 2 

KEY = ’ENTER OUTPUT MACRAN .SCF FILENAME - INCLUDE .SCF : $ 
CALL KEYBD(KEY ,R, I , MACOUT , 3,(12), 0,IR) 

IR = IR + 2 

KEY = ’ENTER OUTPUT MACRAN .SOF FILENAME - INCLUDE . SDF : $ 
CALL KEYBDCKEY, R, I, MACBIN, 3, ( 12), 0,IR) 

C 

C GET THE NUMBER OF POINT IN THE TIME HISTORY 

C 

IR = IR + 2 

KEY = ’ENTER THE NUMBER OF POINTS IN THE TIME HISTORY : $ 
CALL KEYBDCKEY, R, NPTS, A, 2, (4), 0,IR) 

IR = IR + 1 

KEY = ’DO YOU WANT TO TRUNCATE THE NUMBER OF POINTS? $ 

CALL KEYBDCKEY, R,IYN, AYN, 4, (3), 0,IR) 

NEWPT = NPTS 
IF(IYN.EQ.O) GO TO 30 
IR = IR + 1 

KEY =’ ENTER NEW NUMBER OF POINTS : $’ 

CALL KEYBDCKEY , R , NEWPT , A , 2 , ( 4) , 0 , IR) 

30 CONTINUE 

C 

C OPEN UP INPUT COSMOS FILE AND .SCF MACRAN FILE 

C 

OPEN(2,FILE=MACIN) 

C 

C GET INPUT AND OUTPUT CHANNEL ASSIGNMENTS 

C 

IR = IR + 2 

KEY = ’ENTER THE NUMBER OF LOADS TO USE FROM COSMOS : $ 

CALL KEYBDCKEY , R , NFORC , A , 2 , ( 2 ) , 0 , IR ) 

IR = IR + 2 

IF(NFORC.NE.O) GO TO 20 

KEY = ’IS THERE LOAD DATA IN THE COSMOS PLOTFILE? : $ 

CALL KEYBDC KEY ,R, IREWD, AYN , 4 , (3) ,0 , IR) 

IR = IR + 2 
20 CONTINUE 
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c 

C GET CHANNEL DESCRIPTIONS FOR EACH LOAD INPUT 
C 

IF(NFORC.EQ.O) GO TO 8 
CALL VCURXY(O.IR) 

WRITE(*,8020) 

8020 FORMAT( IX, ’ENTER A CHANNEL DESCRIPTION FOR EACH FORCE USED’\) 
CALL VCURXY(O.IR) 

DO 5 IW= 1 , NFORC 
IR = IR + 1 
CALL VCURXY(O.IR) 

KEY = ’ — > $’ 

CALL KEYBD( KEY, R, I , CHANDES( IW),3,(8),0,IR) 

5 CONTINUE 

8 CONTINUE 

CALL VCLEAR 
CALL VCURXY(O.O) 

IR = 2 

KEY=’ ENTER THE NUMBER OF ACCELERATIONS TO USE FROM COSMOS : $’ 
CALL KEYBD( KEY , R , NACC , A , 2 , ( 2 ) , 0 , IR) 

IR = IR + 2 
CALL VCURXY(O.IR) 

WRITE(* , 8030) 

8030 FORMAT( IX, ’ENTER CHANNEL DESCRIPTION FOR EACH ACCELERATION USED’\) 
IP = NFORC + 1 
NTOT = NFORC + NACC 
DO 15 IW=IP,NTOT 
IR = IR + 1 
CALL VCURXY(O.IR) 

KEY = ’ — > $’ 

CALL KEYBD( KEY , R , I , CHANDES (IW),3,(8),0,IR) 

15 CONTINUE 

READ PAST FIRST HEADER IN MACIN, COMPUTE OELTA T AND COMPUTE 


C THE 
C 

NUMBER OF POINTS IN THE FILE 

IR = IR + 2 

CALL VCURXY(O.IR) 

WRITE(*> 7010) MACIN 


7010 

FORMAT( IX, ’READING COSMOS PLOTFILE : 
CALL VCURXY(O.IR) 

REAO( 2,26) ALINE 

’ ,A12\) 

26 

FORMAT (A) 

READ( 2,9010) TIME1 
READ(2,9010) TIME2 


9010 

FORMAT ( IX , E10. 3 ) 

DT = TIME2 - TIME1 
WRITE( * ,7000) 


7000 

FORMAT (80X\) 

CALL VCURXY(O.IR) 
WRITE( * ,7020) DT 


7020 

FORMAT ( IX, ’ FILE HAS A TIME STEP OF : 

’ , 1 PE10. 3 , ’ SECONDS ’\) 


IR = IR + 2 
CALL VCURXY(0 , IR) 

C 

C NOW REWIND FILE AND READ IN DATA 
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REWIND 2 

I F ( (NFORC . EQ. 0) . AND. ( IREWD. EQ. 0) ) GO TO 100 
t IF( NFORC . EQ. 0) GO TO 50 

REA0(2 , 26) ALINE 
WRITER, 7040) 

7040 FORMAT( IX, ’READING IN LOAD DATA FROM COSMOS PLOTFILE ’\) 
IR = IR + 2 
CALL VCURXY(O.IR) 

DO 40 1=1 , NPTS 

READ(2 , 9030) (XIN( IFN , I ) , IFN= 1 , NFORC) 

40 CONTINUE 

GO TO 100 

9030 FORMAT (12X,10(E10.3,1X)) 

9040 FORMAT (12X,3(22X,E10.3,1X)) 

C 

C HERE IF NO FORCES NEEDED BUT MUST BE READ OVER 
C 

50 CONTINUE 

READC2.26) ALINE 
WRITEC* , 7050) 

7050 FORMAT (IX,’ SKIPPING OVER LOAD DATA IN COSMOS PLOTFILE’X) 
IR = IR + 2 
CALL VCURXY(0 , IR) 

DO 60 1=1, NPTS 
READ(2,26) ALINE 
60 CONTINUE 

C 

C HERE TO READ ACCELERATIONS, IF ANY 
C 

100 CONTINUE 

IF(NACC.EQ.O) GO TO 200 
READ( 2 , 26) ALINE 
WRITEC* , 7060) 

7060 FORMAT (IX,’ REAOING IN ACCELERATION DATA’\) 

IR = IR + 2 
CALL VCURXY(O.IR) 

NLINO = (NACC+2)/3 
IP = 1 

103 CONTINUE 

IAC1 = (IP-1) * 3 + 1 + NFORC 
IAC2 = IAC1 + 2 
IF(IAC2.GT.NTOT) IAC2 = NTOT 
DO 105 1=1, NPTS 

REA0(2 , 9040) (XINCIFN , I ) , IFN=IAC1 , IAC2) 

105 CONTINUE 

IFCIP.EQ. NLINO) GO TO 110 
IP = IP + 1 
READ(2 , 26) ALINE 
GO TO 103 
110 CONTINUE 

C 

C HERE WHEN ALL DONE READING IN DATA 
C TIME TO REFORMAT TO MACRAN -.SCF FORMAT 

C 

200 CONTINUE 

WRITEC*, 7070) 

7070 FORMATC IX, ’REFORMATTING TO COSMOS -.SCF FORMAT ’ \ ) 

IR = IR + 2 
CALL VCURXY(O.IR) 

NPTS = NEWPT 
CALL OUT 1 
C 
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C HERE TO WRITE OUT REST OF MACRO 

C 

OPEN (6 , FI LE= ’ NASAMAC . BAT ’ , STATUS= ’ NEW’ ) 

WRITE(6, 19) 

19 FORMAT (’ INTRFS.EXE’) 

WRITE(6 , 1911) MACOUT 

1911 FORMAT (’MAC42 3 /R 10000 < ’ ,A12, IX, ’ >MACOUT.OUT’ ) 
WRITE(6, 1912) 

1912 FORMAT (’DEL MACOUT. OUT’) 

WRITEC6, 1913) 

1913 FORMAT (’NASADAM’) 

CLOSEC6) 

HERE WHEN DONE 

STOP 
END 
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MODAL STRAIN ENERGY DISTRIBUTION METHOD (MSEDM) 
SOURCE CODE LISTINGS 
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INOFLOATCALLS 

IMPLICIT INTEGER*2 (I-N) 

DIMENSION DAM ( 20) , DAMCRT (20) 

CHARACTERS KEY 
CHARACTER *60 DAMTTL 
CHARACTER*20 DAMFIL 
CHARACTER* 1 IEDCR 
C 

C CLEAR SCREEN, ETC 
C 

CALL VCLEAR 
CALL VCURXY(O.O) 

KEY = ’ DAMIO - PROGRAM TO CREATE/EDIT DAMAGE PARAMETERS$ 
CALL PRMPT( KEY ,1.00) 

IR = 6 
IR = IR + 1 

KEY = ’DO YOU WANT TO (C)REATE OR (E)DIT A DAMAGE FILE? : $ 
CALL KEYBD( KEY , R, I , IEDCR ,3,(1),0,IR) 

C 

C HERE FOR FILE CREATION 
C 

CALL VCLEAR 
IR = 0 

CALL VCURXY (0 , IR) 

KEY = ’DAMAGE TITLE : $’ 

CALL KEYBD(KEY,R, I .DAMTTL, 3, (60) ,0, IR) 

IR = IR + 1 

KEY = ’NUMBER OF MODES USED : $’ 

CALL KEYBD(KEY , R , NMODES , A , 2 , ( 2) , 0 , IR) 

IR = IR + 2 
CALL VCURXY(O.IR) 

WRITE(* , 9010) 

9010 FORMAT (IX, ’ MODE ’ , 5X , ’ D(W) ’ , 1 0X, ’ Deri t(W) ’\) 

IR = IR + 1 
CALL VCURXY(0, IR) 

WRITE(* , 9002) 

9002 FORMAT (’ N0.’\) 

IR = IR + 1 
C 

C NOW READ IN DATA 
C 

1ST = IR 
I END = IR + 9 
DO 20 IM=1, NMODES 
CALL VCURXY(O.IR) 

IF(IR.EQ.IST) WRITE(* , 9200 ) IM 
IF(IR.NE.IST) WRITE( * , 920 1 ) IM 
KEY = ’ $ ’ 

CALL KEYBD(KEY,DAM(IM),I,A, 1 ,(10),9,IR) 

CALL KEYBD(KEY , DAMCRT ( IM) , I , A , 1 , (10),23,IR) 

IR = IR + 1 

IF(IR.LE.IEND) GO TO 20 
DO 15 IER= 1ST , IEND 

CALL VCURXY(0 , IER) 

WRITE( * , 9000 ) 

15 CONTINUE 

IR = 1ST - 
20 CONTINUE 
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C 

C GET OUTPUT FILENAME AND WRITE IT OUT 

C 

KEY = ’ENTER OUTPUT DAMAGE FILENAME : $ 

CALL KEYBD(KEY,R , I , DAMFIL ,3,(20),0,IR) 

OPEN( 5 , FILE=DAMFIL , STATUS= ’ NEW’ ) 

WRITE( 5 , 8000) DAMFIL , DAMTTL , NMODES 
WRITEC 5 ,8010) 

DO 30 1=1, NMODES 

WRITEC 5 , 8020 ) I , DAM( I ) , DAMCRT (I ) 

30 CONTINUE 

9201 FORMAT (IX, 1X,I3\) 

9200 FORMAT (IX, I 3\) 

9000 FORMAT (40X\) 

8000 FORMAT (’DAMAGE FILENAME : A20/A60/ ’ NUMBER OF MODES 

8010 FORMAT ( ’ MODE ’ ,5X, ’D(W)’ , 10X, ’Dcrit(W)’ ) 

8020 FORMAT ( 1 X , 12 , 5X , F7 . 3 , 7X , F7 . 3 ) 

STOP 

END 
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SUBROUTINE MODSSUM 
IMPLICIT INTEGER*2 (I-N) 
DIMENSION ESC 100 , 20) 
CHARACTERS KEY 
CHARACTERS TITLE 

C 

C GET CURSOR 

Q 

CALL GCURXY ( I C , I R ) 

IR = IR + 1 


C OPEN UP FILE 

U READC5) TITLE, NELM, NMODES 

DO 10 1=1, NMODES 

READ(5) ( ES ( J , I ) , U~ 1 i NELM ) 

10 KEY = ’ENTER WHICH MODE TO SUMMARIZE : $’ 

CALL KEYBDC KEY , R , IMODE , A , 2 , ( 3) , 0 , IR) 

IR = IR + 1 
C 

C OPEN UP MODSUM.DD FILE 

^ OPEN(6 , FILE= ’ MODSUM. DO ’ , STATUS= ’ NEW’ ) 

WRITE(6 , 8000) TITLE, NMODES, NELM 

8000 FORMAT (A40, 

1 /’NUMBER OF MODES : ’,12, 

2 /’NUMBER OF ELEMENTS : ’,13) 

WRITE(6,8010) IMODE 

8010 FORMAT (//20X, ’SUMMARY FOR MODE NUMBER : ,12) 

WRITE(6 , 8020) M . niTU , 

8020 FORMAT (//’ ELEMENT ’,5X,’ MODAL STRAIN , 

1 /’ NO. ’ , 5X , ’ ENERGY’) 

DO 20 1=1 ,NELM 

WRITE(6,8030) I , ES( I , IMODE) 

20 CONTINUE 

8030 FORMAT ( 3X , 1 3 , 9X , 1 PE 1 0 . 3 ) 

STOP 

END 
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C 

C PROGRAM MOOSTRN - COMPUTE MODAL STRAIN ENERGY MODE BY MODE 

C 

IMPLICIT INTEGERS (I-N) 

DIMENSION XTRN(200),XAVG(200) 

CHARACTER*80 KEY 
CHARACTER*40 TITLE 
COMMON/XDAT/XTRN , XAVG 


HObsTfc^. F0<^ 


c 

C CLEAR SCREEN, DISPLAY PROGRAM 
C 

KEY = ’MOOSTRN - PROGRAM TO COMPUTE MODAL STRAIN ENERGY $’ 
CALL VCLEAR 
CALL VCURXY(O.O) 

CALL PRMPT (KEY ,1.00) 

CALL GCURXY ( IC, IR) 

IR = IR + 1 


C 

C GET USER INPUTS 
C 

KEY = ’ENTER DESCRIPTION OF RUN : $’ 

CALL KEYBD(KEY, R , I ,TITLE, 3, (40) ,0, IR) 

IR = IR + 1 

KEY = ’ENTER THE NUMBER OF ELEMENTS TO USE : $’ 
CALL KEYBD(KEY,R,NELM,A,2, (2) ,0,IR) 

IR = IR + 1 

KEY = ’ENTER THE NUMBER OF MODES TO USE : $’ 
CALL KEYBD( KEY , R , NMODES , A , 2 , ( 2 ) , 0 , IR) 

IR = IR + 1 
CALL VCURXY(O.IR) 

C 

C OPEN UP BINARY FILE CONTAINING DATA 
C 

0PEN(6 , FILE= ’ MS. DD’ , STATUS= ’NEW’ , 

1 FORM=’ UNFORMATTED’) 

WRITE(6) TITLE, NELM, NMODES 

C 

C COMPUTE STRAIN ENERGY FOR 100 PTS FROM 0 TO 1 
C 

XLM = 1.0 / NELM 
NELMP1 = NELM + 1 
XPI = 4.0 * ATAN(I.O) 

C 

C LOOP OVER NUMBER OF MODES 
C 


DO 100 IM0DE=1 .NMODES 
XMX =0.0 
XCOFF = 1.0 
DO 10 1=1, NELMP1 

XL = FLOAT ( 1-1 ) * XLM 
XS = FLOAT (IMODE) * XPI 
XSN = SIN(2*XS*XL) 

XTRN(I) = XCOFF * (XL/2 - XSN/(4.0 * XS)) 
10 CONTINUE 

C 

C COMPUTE ELEMENT STRAIN ENERGY 
C 
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4 DO 15 1= 1 , NELM 

XAVG(I) = (XTRNU + 1) - XTRN( I ) )/0. 1 
IF(ABS(XAVG(I)).GT.XMX) XMX = XAVG(I) 
15 CONTINUE 

SCALE TO UNITY 

DO 20 1=1, NELM 

XAVG(I) = XAVG(I) / XMX 
20 CONTINUE 

NOW WRITE IT OUT 

WRITEC6) (XAVG(I), 1=1 , NELM) 

100 CONTINUE 

STOP 
END 
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i $NOFLOATCALLS 

IMPLICIT INTEGER*2 (I-N) 

DIMENSION ES( 1 00 , 20) , DAM ( 20) , DAMCRT (20) 

CHARACTER*80 KEY 
CHARACTER*60 DAMTTL 
CHARACTER*40 TITLE 
CHARACTER*20 FLNME , DAMFIL 
CHARACTER*3 P0AM( 100) , ISM 
CHARACTER*27 FMT 
CHARACTER*2 FM 

CHARACTER* 1 AES( 100 , 20 ) , FF, FMTM( 27 ) , FMM(2 ) 

DIMENSION ESMX(20) , LESMX(20) 

COMMON/DAT/ES , AES , ESMX , LESMX , DAM , DAMCRT , PDAM 
EQUIVALENCE ( FMT , FMTM( 1 ) ) , ( FM , FMM( 1 ) ) 

DEFINE FORMFEED 

FF = CHAR(12) 

FMT = ’ ( 7X , 12 , 5X , 1(A1 ,3X) , 10X.A3)’ 

FM = ’ T 

DEISPLAY PROGRAM TO USER 

CALL VCLEAR 
CALL VCURXY(O.O) 

KEY = ’MODSUM - PROGRAM TO SUMMARIZE STRAIN ENERGY DATA$ ’ 

CALL PRMPT( KEY ,1.00) 

CALL GCURXY ( I C , I R ) 

IR = IR + 1 

GET FILENAME CONTAINING DATA 

KEY = ’ENTER FILENAME CONTAINING MOOAL STRAIN ENERGIES : $’ 
CALL KEYB0( KEY , R, I , FLNME ,3,(20),0,IR) 

OPEN( 5 , FILE=FLNME , FORM= ’ UNFORMATTED’ ) 

IR = IR + 1 

SEE IF USER WISHES TO SUMMARIZE MODAL STRIN ENERGIES 
FOR ALL ELEMENTS 

KEY= ’ DO YOU WANT TO SUMMARIZE STRAIN ENERGIES FOR A MODE? : $’ 
CALL KEYBD( KEY , R, I , ISM ,4,(3),0,IR) 

IR = IR + 1 

IF(I.EQ.I) CALL MODSSUM 

GET FILENAME CONTAINING DAMAGE INFO 

KEY = ’ENTER FILENAME CONTAINING DAMAGE DATA : $’ 

CALL KEYBD(KEY,R, I , DAMFIL , 3,(20),0,IR) 

IR = IR + 1 
OPEN (7 , FI LE= DAMFIL) 

GET DAMAGE FRACTION 

KEY = ’ENTER DAMAGE FRACTION ( <= 1 ) : $’ 

CALL KEYBD(KEY,FR, I, A, 1 , (10) ,0, IR) 

C 

C READ IN DATA 
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READ ( 5 ) TITLE, NELM,NMODES 
WRITECFM, 9111) NMODES 
FMTM(II) = FMM( 1 ) 

FMTM(12) = FHM( 2 ) 

9111 FORHAT(I2) 

READ(7 ,7010) DAMTTL 
7010 FORMAT (/A60/ 1 X/1X) 

DO 10 1=1, NMODES 

READ(5) ( ES( J , I ) , J= 1 , NELM) 

READ(7 ,7020) DAM( I ) , DAMCRT ( I ) 

10 CONTINUE 

7020 FORMAT (8X , F7 . 3 , 7X , F7 . 3) 

C 

C NOW DETERMINE MAX STRAIN ENERGIES FOR EACH MODE 
C 

DO 30 1=1, NMODES 
XMX =0.0 
DO 20 J= 1 , NELM 

IF(ES(J,I). LE. XMX) GO TO 20 
XMX = ESC J , I ) 

LMX = J 


20 

CONTINUE 



ESMX(I) = 

XMX 


LESMX(I) 

= LMX 

30 

CONTINUE 



C 

C INITIATE AES 
C 

DO 35 1=1, NMODES 
DO 34 J= 1 , NELM 

AES( J , I ) = ’ ’ 

34 CONTINUE 

35 CONTINUE 
C 

C NOW SORT IT ALL OUT 
C 

C BEGIN BY FINDING ELEMENTS WITH LOW STRAIN ENERGIES 
C 

DO 100 1=1, NMODES 

I F ( DAM ( I ) . GT . DAMCRT ( I ) ) GO TO 50 
ESMXL = FR * ESMX(I) 

DO 40 J= 1 , NELM 

IF(ES(J,I).LE. ESMXL ) AES(J,I) = ’x’ 
40 CONTINUE 

GO TO 100 
C 

C HERE FOR ELEMENTS WITH HIGH STRAIN ENERGIES 
C 

50 CONTINUE 

DO 60 J=1 , NELM 
AES( J , I ) = ’x’ 

60 CONTINUE 

100 CONTINUE 

C 

C NOW DETERMINE IF ANY ONE ELEMENT MIGHT HAVE DAMAGE 
C 

DO 110 IL= 1 , NELM 
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PDAM(IL) = ’YES’ 

DO 120 IM= 1 , NMODES 

IF(AES( IL , IM) . EQ. ’ ’) PDAM(IL) = ’ NO’ 

120 CONTINUE 
110 CONTINUE 

C 

C NOW SUMMARIZE DATA 

C 

OPEN ( 6 , FILE= ’ MODSUM. DD ’ ,STATUS= ’NEW’ ) 

NPG = (NELM + 39 )/40 + 1 
I PAGE = 1 

WRITE( 6 , 9000 ) TITLE , IPAGE , NPG , DAMFIL , DAMTTL , NMODES , NELM , FR 
9000 FORMAT ( //4X, A40, 16X, ’Page ’,11,’ of \I1, 

1 /4X, ’DAMAGE FILE USED : \A20 

1 /4X, ’DAMAGE TITLE : ’ ,A60 

1 /4X, ’NUMBER OF MODES : ’,12, 

2 /4X, ’NUMBER OF ELEMENTS : ’,12, 

3 /4X, 'DAMAGE FRACTION : ’,F6.4) 

WRITE(6 , 9010) 

9010 FORMAT (/4X, ’MODE’ ,3X, ’MAXIMUM ELEMENT 3X, ’ ELEMENT ’ , 

1 /4X , ’ NO. ’ , 3X , ’ STRAIN ENERGY ’,3X,’ NO.’,/) 

DO 200 1=1, NMODES 
WRITE(6 , 9020 ) I , ESMX(I) , LESMX(I) 

200 CONTINUE 

9020 FORMAT ( 5X , 1 2 , 6X , 1 PE 1 0 . 3 , 8X , 1 2 ) 

WRITE(6 , 902 1 ) FF 

9021 FORMAT (A1) 

IPAGE = 2 

WRITE ( 6 , 9000 ) TITLE, IPAGE , NPG , DAMFIL , DAMTTL , NMODES , NELM , FR 
WRITE(6 , 9035 ) 

9035 FORMAT (//20X, ’ELEMENT DAMAGE SUMMARY’) 

WRITE(6,9040) 

9040 F0RMATC/4X, ’ELEMENT’ ,27X, ’MODE’ ,27X, ’POSSIBLE’ , 

1 /4X , ’ NO. ’ , 27X , ’ NO. ’ ,28X, ’DAMAGE’ ) 

WRITE(6 , 9050 ) ( I , 1= 1 , NMODES) 

9050 FORMAT(13X, 15(12, 2X)) 

DO 70 J= 1 , NELM 

WRITE( 6 , FMT ) J , (AES(J , I ), 1= 1 , NMODES) , PDAM(J ) 

IF(MOD( J , 40) . NE . 0) GO TO 70 
WRITE(6 , 9060) 

IPAGE = IPAGE + 1 
WRITE(6 , 902 1 ) FF 

WRITE(6 , 9000) TITLE , IPAGE , NPG , DAMFIL , DAMTTL .NMODES , NELM , FR 
WRITE(6 , 9035 ) 

WRITE( 6 , 9040 ) 

WRITE(6 , 9050 ) ( I , 1= 1 , NMODES) 

70 CONTINUE 

WRITE(6 , 9060) 

9060 FORMAT (// 1 5X, ’x Indicates possible modal element damage’) 

C 

C ALL DONE 
C 

STOP 

END 
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FOREWORD, 


The MAC/RAN IV System is a proven program for the analysis of time 
series data. Originally developed by Measurement Analysis 
Corporation of Los Angeles, MAC/RAN was acquired by University 
Software Systems in 1970. The program has been continuously refined 
and improved since that time and has served data processors and 
analysts in a wide variety of applications. 

MAC/RAN IV differs from MAC/RAN III in a number of ways. Once an 
industry standard in the mainframe computer world of the 1970s, the 
advent of supermini, mini, and microcomputers led to the complete 
restructuring of MAC/RAN III. MAC / RAN IV can now take full advantage 
of the greater high speed memory (if available) of modern computers, 
virtual memory operating systems, and fast mass storage devices. A 
new, more compact data file structure, along with its direct access 
allows MAC/RAN to run efficiently on desktop microcomputers as well as 
the largest of mainframes. The Data File Manager controllers 
(replacing the EXEC "cards" of MAC/RAN III) reflect this change in 
data access and structure. 

The first three chapters of this manual describe the terms and 
concepts used in MAC/RAN IV along with the structure of the program. 
Chapter 4 describes the Data File Manager. Chapters 5 through 15 
describe the MAC/RAN IV processors. It is recommended that the user 
first carefully read Chapters 1 through 4 and study the examples (and 
possibly the Standard Test Case of Appendix B) to get a feel for the 
processing sequences. Chapters 5 and up may be referenced as needed. 

The PREP (Chapter 5) and PSD (Chapter 8) Processors are particularly 
powerful and offer many attractive features. PREP contains the 
digital filters and PSD does spectral averaging (as opposed to simple 
frequency transformation in FOUR) along with coherency, cross 
spectra, and transfer function computation. Each computational 
processor has a section that describes the method or procedures that 
it uses . 

MAC/RAN generally operates in a batch mode and is, in fact, designed to 
process large streams and many channels (up to 25 at one time) of data 
efficiently and with minimum user interference. However, if 
desired, MAC/RAN can now operate on a controller by controller basis, 
directly from the computer console. The main differences between 
this mode of operation and a truly interactive one. is that the user 
must have a grasp of the control language as MAC/RAN does not supply 
menus or help lists; the reference manual provides the information. 
The other difference from truly interactive operation is that the 
plots generated are not immediately displayed, but are written to a 
file for later postprocessing or output, depending on equipment and 
installation. Thus MAC/RAN IV remains essentially hardware 
independent during its operation. 


iii 
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CHAPTER 1 


INTRODUCTION 


A. BACKGROUND 

1 . Tiae Series Analysis 

A series of data values that are compiled with respect to their 
chronological occurences can be considered a tiae series. This 
■eans that a data value or perhaps an enseable of values are 
associated with one instant in tiae that in aost cases precedes 
and succeedes another pairing or enseable of values. A siaple 
exaaple sight be recording the value of several stocks on a daily 
basis . 

Tiae series analysis Beans applying analytic aethods to better 
understand tiae histories of values or to reduce thea into a fora 
wheTe conclusions can aore readily be aade froa the data. Soae 
typical exaaples of tiae series analysis aight be plotting the 
data values with respect to their tiae occurrences, liaiting the 
rates at which these values change (saoothing or filtering), or 
correlating one set of data points with anothar. 

More often, one is concerned with a transforaation of tiae series 
data into another doaain such that aore involved investigation 
can be facilitated. The aost coaaon of these transf oraat ions is 
the Fourier transfora which in soae sense calculates the 
frequency content of a tiae history of values. This 
transforaation is a very powerful tool and is used extensively 
for aany purposes, including power estiaation, spectral display 
and analysis. certain types of correlation, and even 
aatheaatical integration and differentiation, aaong aany 
others. Thus, tiae series analysis is a very broad tera that 
includes aany frequency doaain techniques to deal with a wide 
spectrua of phenoaena. These techniques are also applicable to 
situations that do not involve a tiae basis, since an analogy can 
be aade with a tiae series. 

2 . Solving Tiae Series Probleas 

A typical way to solve isolated tiae series probleas is to 
develop individual, non- i ntegrated prograas. This usually 
requires soae research into the techniques eaployed, a choice of 
algorlthas, and then the tiae consuaing tasks of prograaaing and 
debugging on a coaputer. With an arbitrary collection of 
prograas to do the different kinds of signal and data analyses 
required, probleas of coapat i bi 1 i ty are likely. For exaaple, 
output froa one prograa aay not have the proper foraat for input 
to another prograa. Thus, prograas and/or data nay have to be 
partially rewritten nearly every tiae they are used. Frequent 
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program and data changes increase the probability for human 
error. Corrections are costly and time-consuming; more 
significantly, undetected errors may cause distorted or 
incorrect results. 

3 . mac/RAN as a Time Serie s_Tooi 

MAC/RAN incorporates sophisticated algorithms and data 
preparation techniques in its processors. This minimizes the 
signal processing expertise and computer coding necessary to 
analyze data correctly. All data entered into MAC/RAN remains 
in a standard data format that allows continuous processing 
through various modules or selective output at any time during a 
program run. In addition MAC/RAN processes multiple channels 
simultaneously for greater throughput. 

NAC/RAN allows any one, several, or all of its state-of-the-art 
processing procedures to be specified within the same program 
run. The MAC/RAN driver and management software handle the 
delegation of tasks and data to the necessary modules by way of a 
command language. Thus, a command file can be easily edited to 
alter existing sequences of procedures, parameters, or input 
data. This avoids the rewriting of special software that is 
required to accommodate any adjustments needed to obtain the 
desired results. 

MAC/RAN is not a programming language nor a set of subroutines 
from which one constructs his own data reduction system; 
instead, it is a self-contained data reduction system that 
provides flexibility in the choice of any specific processing 
procedures required at any given time. 


B. SYSTEM FEATURES 
1 • F lex ibi 1 i ty 

As a package of state-of-the-art processing procedures. MAC/RAN 
performs a wide variety of basic computations to properly reduce 
time series data. Any one, several, or all of the basic 
processors (plus the Data File Manager) may be specified within 
the same program to solve a given problem. Figure 1-1 
illustrates the MAC/RAN system. 

Processors include: 

- Data File Manager. (FILIN, FILOUT, MERGE) 

- Data Preparation Processor (PREP) 

- Fourier Analysis Processor (POUR) 

- Power Spectral Density Processor (PSD) 
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- Plugboard Simulation Processor (PLUG) 

- Display Processor (PLOT) 

- Input/Output Port (CTAPE) 

- Convolution and Filtering Processor (CONFIL) 

- Time and Frequency Analysis Processor (SPEC) 

- Amplitude Statistics Processor (AMSTAT) 

- Shock Spectrum Analysis Processor (SHOCKS) 

- Linear Systems Analysis Processor (TRANS) 



Figure 1-1. MAC/RAN System 
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2 . S lap 1 1 f led Programing Requiremen ts 

This modular construction of MAC/RAN provides a simplified 
approach to processing tiae series data. The Data File Manager 
controls the flow of data from processor to processor in order to 
compute the functions specified in any given run. This feature 
significantly reduces the aaount of control information the user 
must supply to link each specified processor. 

3 . integrated Sy s tem Capability 

Through MAC/RAN's package of general programs, the complete 
processing for a given data set of time histories can be obtained 
with a single computer run; in contrast, a non- i ntegrated system 
would extend the processing over a long period of time with 
individual computer runs required for each of the individual 
functions. However, the independence of the processors, even 
though they can be combined into a single package, allows any 
given total processor operation to be repeated any number of 
times in a single run. This flexibility in control of processing 
can be implemented by the control language parameters. 

4 . Dat a Format Compatibility 

Throughout the system, every intermediate data format is the 
same; thus, input to any given processor is in the same format. 
Consequently, after data is initially fed into the MAC/RAN 
System, it need never be manually handled again. This greatly 
simplifies bookkeeping details and increases the efficiency of 
the overall processing procedure. 

5 • Displa y Capability 

A special processor is devoted to providing printed and/or 
plotted outputs of any data or any intermediate functions 
stored in the standard MAC/RAN format. Inherent within the 
display processor is the ability to select logarithmic or 
decimal scales, plus the basic linear scale and scale increment 
selections in numerically simple values. 


C. SYSTEM CONFIGURATION 

The MAC/RAN General Process Chart, shown in Figure 1 -2 , i 1 1 ustr ates 
the flow of data through the system, and the interrelationship between 
the various stages of processing. These stages are described by 
symbols representing software, hardware and people. 
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1 . I nout 

Input to MAC/RAN consists of a control file and data to be 
analyzed for the desired processors. 

2 . Contro l 

Control of processing is aaintained by the Systea Manager, 
invisible to the user, through the paraaeters in the control 
file. 

3 . Outp ut 

Generally, a coaputational processor prints its operational 
inforaation and the channel statistics of the output data. This 
processed data is written out as a Standard Data Pile (SDP) . The 
Display Processor (PLOT, not considered coaputational) aay be 
used at any tiae to print or plot data froa the SDP of any MAC/RAN 
processor. Soae coaputational processors optionally print and 
plot their own output during execution. 
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Figure 1-2. MAC/RAN General Process Chart 
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CHAPTER 2 - SYSTE M OPERATION 


A. GENERAL 

MAC/RAN consists of a Systea Executive, a Data File Manager and 
various processors that do the tiae series analysis. In addition, 
there are two special processors that handle I/O and coaaunicate 
MAC /RAN results to the user, respectively. User Inpu ts to M AC/RAN 
are a Standard^ Control _F_1 le_ ( S_C F. Land (in aost cases) one or aore data 
fi les . These data files can contain aul tiple channels of (Tat iT that 
are to b^~lTFbcYs~sed s iauTtaneous 1 y , if desired. 


The systea operation follows the user coaaands placed in the input 
control file. The Systea Executive recognizes these user specified 
coaaands and transfers control to the data file aanager or one of the 
processors . 

The MAC/RAN Executive takes care of systea initialization, control 
file inforaation, calls the required processors, and constantly 
updates the channel statistics during processing. During the course 
of the processors' coaputat ions , the Executive also decides when to 
write the output data to the appropriate disk file. 

The Data File Manager can coabine input data files, select or coablne 
channels, and edit the data to be processed. Editing features are 
start and stop tiae specification, leading or trailing zero padding, 
specification of nuaber of points desired. and deciaation. 
Interpolation is possible using the PREP Processor. 


B. STANDARD CONTROL FILE (SCF) FORMAT 

A Standard Control File consists of various control records, or rows 
of control inforaation, froa which the Systea Executive drives 
MAC/RAN. These control records Bust confora to a standard foraat of 
seven fields per record. 

The first three fields are alphanuaeric and generally contain 
processor naaes and processor functions (controllers) as well as 
annotation inforaation. Field 2 is used only for continuation with 
certain processors. Fields 4 through 7 are nuaeric (except when the 
JBHEAD controller is used). Nuabers in these fields can be placed 
anywhere if a deciaal point is supplied, otherwise they aust be right 
justified. Generally, blanks in nuaeric fields are interpreted as 
zeros. A coaplete description of the specific control record foraats 
are given in the following chapters. 

The control field foraat is shown below. The arrows aark the 
beginning of the field that follows. 
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Control Field Poraat: 

I 7a 24 38 48 80 7| 

i /; II n i f± i /$ * /g i IZ 

Record 1 
Record 2 


C. MAC/RAN STANDARD DATA PILE (SDP) FORMAT 

All data files generated by the systea are written in a standard binary 
foraat. The first record of a SDP stores channel naaes and statistics. 
The subsequent records store channel data. The channel statistics 
consist of the aean value, variance, aaxiaua value, ainiaua value, 
nuaber of channel data points, and saapling rate. 

All MAC/RAN processors generate channel statistics when they output a 
SDP file. The channel (or function) naaes are defined by the user and 
they provide annotation for plots. If no naaes are entered, MAC/RAli 
supplies then as DATA 1 , DATA2 , etc. If the naaes are properly chosen, 
then the plots will be easily identifiable. 


D. DEFINITION OP TERMS EMPLOYED 
1 . General 


Controller - A MAC/RAN control word or a part of a control 

record that specifies control inforaation to 
the systea; e.g.. "PLUG" or "FILTLOW" (in the 
PREP Processor). It aust be capitalized . 

Control Field - One of seven specific areas of a control record 

where MAC/RAN control paraaeters are placed. 


Control Pile 


A disk file that contains the control records 
to drive MAC/RAN in its specific tasks. 


Control Group - Adjacent control records with the sane 

controller in a MAC/RAN control file. 


Control Record - A record in the MAC/RAN control file that 

contains a controller and optionally, other 
control paraaeters. 


DFM 


The Data Pile Manager. 
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DTSA 


Praae 


Pune t ion 


"Digital Tiae Series Analysis” by Otnes A 
Enochson ( W i 1 ey- I n ter sc ience , 1972). 

A data group coaprlsing one eleaent of the 
independent channel and the corresponding 
function (dependent channel) values that 
occur s iau 1 1 aneous 1 y . e.g. a tiae slice froa a 
data f il_e that consists of lO'channeia - ; One 
frame would contain 10 elements.-- tiae and nine 
function values. 

The data that resides in a dependent channel. 


SCP 


A MAC/RAN Standard Control Pile. 


SOP - A MAC/RAN Standard Data Pile, containing 

fraaes of data and channel statistics. Each 
channel (function) in the SDP has the saae 
length. 


2 . Matheaat leal 


A tiae series x, will be assuaed to be saapled at equal spaced 
I nterval s seoar a t ed in tiae bv an amoun t — A-tr 30 that the 
continuous fora of x,. x(t), is related to it by 

x ( t ) - x( iAt ) - x, (2.1) 

Without loss of generality, i will be assuaed to have the range 
0,1,.., N- 1 for a total of N values. 

The series y ( will frequently be used to represent either a 
second series, or the original series x, after soae operation 
such as digital filtering. 

Pourier transforas will be denoted with capital letters. Por - 
exaaple, the Pourier transfora of x(t) will be written as X(f) 
w here f is frequency in Hz . The Pourier transform of x, will be 
written X k where 

f k - kAf (2.2) 

and Af is the frequency interval in Hz. Usually. 

1 

f k » (2.3) 

NAt 

The saapling rate f s is usually input to the program rather than 
At, but their relation is quite simple 
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1 

At - (2.4) 

f 3 

The aaxiaui frequency. f N , soaetiaes referred to as the Nyquist 
frequency is 


f 3 

f M - « — 

2 2 At 


(2.5) 
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PLOT 


CHAPTER 12 - DISPLAY PROCESSOR (PLOT) 


A. FUNCTION 

The Display Processor provides a convenient leans for listing and/or 
plotting output froa any data file in the NAC/RAN SDF foraat. Any 
function or coabination of functions fro* a file coaputed with the 
NAC/RAN Systea and written on disk in the SDF foraat aay be saved and 
later printed or plotted by the Display Processor. 


B. INPUT/OUTPUT SPECIFICATIONS 

1 . Input 

a . Control Inputs to the Display Processor include: 

o Printer and/or plot selection 
o Multiple plot selection 

o Plot widths and lengths 

o Plot data intervals 

o Plot scale selection — linear, log, or db 

b . Control Inputs to Data File Manager ( DFM ) . The FILIN 
controller (and optional FILIN CCHAN Controller) is 
required to describe input disk files. 

c. Data Inp ut. Input aay consist of one or lore channels of 
time series data in MAC/RAN SDF foraat. 

2 . Out put 

Any desired listing and/or plots with one to 24 curves per plot 
aay be specified. As indicated above, SDPs that have been 
produced by any coaputat i onal processor in standard foraat aay 
be reprinted or replotted at any tiae. 


C. CONTROL RECORDS 

Control records for the Display Processor aay be input in any order 
with the exception of aultiplot continuation records which aust 
follow the aultiplot definition record (see C.l.g below). All 
control records are optional. 

1 . Forma t 

The basic foraat of each record is shown in the following figure. 
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1 7« 

i ft 


PLOT PLOT 

PLOT PRINT 

PLOT PRINTIND 

PLOT INDE II 

PLOT LOG T 

PLOT SCALE VV T 

PLOT NAMEXXXXXXXXXXXX 

PLOT MULTa II NN T 

PLOT CMULT 

PLOT CNANEXXXXXXXXXXXX 


XCYC 

XO 


YCYC 

XMAX 

V2 

YO 

V 3 

YMAX 

V 4 

XNAX 

V2 

YO 

V 3 

YMAX 

V 4 


Plot/Print Record 
Field 1 - PLOT 

Field 3 - Print and single plot output options 

Blank 4 Print and plot 

PLOT •» Plot only - (left adjusted) 

PRINT 4 Print only - (left adjusted) 

Individual Print Selection Record 

Field 1 - PLOT 

Field 3 - PRINTIND 

If print is requested through Record (a), a simultaneous 
printout of all input data channels is nominally produced. 
If the Individual Print Selection Record is input, however, 
a separated printout of the independent variable along with 
each dependent variable is produced rather than the 
simultaneous printout. 

Independent Variable Selection Record 
Field 1 - PLOT 

Field 3 - INDE - (left adjusted) 

Independent channel number (numerical 1-25) - 

(columns 16-17. right adjusted) 

The independent data channel for printing and for the 
individual plots is nominally Channel 1. This may be 
overridden by specifying the Independent variable through 
this control ler . 
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PLOT 


i 


d. 


Plot Type and Log Plot Cycle Size Selec t ion R eco r d 


Field 1 - PLOT 

Field 3 - LOG - (left adjusted) 


Nominal pl^ot type (numerical 1-8) 

1 + Linear x. linear y 

2 + Linear x. log y 

3 •* Log x. linear y 

4 ■» Log x. log y 

5 ■> Linear x. 10 log t Q(y) ' 

6 ♦ Linear x. 20 log ( g(y) 

7 *► Log x. 10 1 og j q ( y ) 

8 + Log x, 20 log|o(y) 


- (column 23) 


db 


The nominal plot type (type 1 unless specified by 
this record) is used for all plots for which no 
specific plot type is requested (see the SCALE 
record and MULT record below). 

Field 4 - Log cycle size in inches for x-axls (floating 
point). If not specified, 2.0 is used. 

Field 5 - Log cycle size in inches for y-axis. If not 
specified, 2.0 is used. 


e • Scale Selection Record 
Field 1 - PLOT 

Field 3 - SCALE - (left adjusted) 

Dependent channel number, dependent variables 
only (numerical 2-25). (Columns 16-17, right 
adjusted.) Defines variable number for which 
the Scale Selection Record parameters apply. 
If this field is blank, the record parameters are 
applied to all individual plots. 

Plot Type - (numerical 1-8) - (column 23). 

Specified plot type is used for plots implied by 
above channel number field. If blank, nominal 
plot type is used (see C.l.d above). 

Field 4 - Origin for x-axis (floating point) 

Field 5 - Maximum value for x-axis (floating point) 

If Fields 4 and 5 are not specified, processor 
selects x-axis scales. 
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Field 6 - Origin for y-axis (floating point) 

Field 7 - Naxleua value for y-axis (floating point) 

If Fields 6 and 7 are not specified, processor 
selects y-axis scales. 

f . Pl ot Axis Naee R e cord 

Field 1 - PLOT 

Field 3 - NAME (left adjusted) 

12-character naee field (columns 12-23) 

Specifies a 12-character label for plot axes for 
channel numbers given in Fields 4 through 7. If 
not specified for a channel, the name (ID) field 
on the data file is used. 

Fields - Channel numbers to mhich specified name field 
4-7 applies. 

* • Multiple Plot Definition Records 

The Multiple Plot Definition Records define a single plot 
containing multiple dependent variables plotted against a 
single Independent variable on the same set of axes. Only 
one multiplot may be requested in one execution of the 
Display Processor. 

Field 1 - PLOT 

Field 3 - MULT (left adjusted) 

Independent Channel Number - (numerical 1-25) - 
(columns 16-17, right adjusted). Specifies 
independent variable for multiplot. If 
omitted, the nominal independent variable is 
used (see C.l.c). 

Number of dependent channels to be plotted - 
(numerical 1-24) - (columns 20-21. right 

adjusted). If blank, all data channels except 
the independent channel are plotted. 

Plot Type - (numerical 1-8) - (column 23). If 
omitted, the nominal plot type is used (see 
C. 1 .d) . 

Field 4 - Origin for x-axls (floating point) 
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Field 5 

Field 6 
Field 7 


- Maximum value for x-axis (floating point) 

If Fields 4 and 5 are not specified, processor 
selects x-axis scales. 

- Origin for y-axis (floating point) 

- Maximum value for y-axis (floating point) 

if Fields 6 and 7 are not specified, processor 
selects y-axis scales. 

1 ) Multiplot Continuation Record 

Field 1 - PLOT 

Field 2 - C 

Field 3 - MULT (left adjusted) 

Fields - Dependent data channels to be 
4-7 plotted - (numerical 1-25) 

columns 34-35 , 46-47, 58-59, 70- 
7 *)- If specific dependent 
data channels are not specified, 
the first n channels (not 
including the Independent chan- 
nel) are used where n is 
specified by the Multiple Plot 
Definition record above. If 
fewer than n channels are 
explicitly requested, the re- 
mainder are automatically se- 
lected starting from the first 
variable (omitting the indepen- 
dent channel and avoiding 
duplication). As many Multi- 
plot Continuation Records as 
necessary may be input. 

2 ) Multiplot Name Record 

Field 1 - PLOT 

Field 2 - C 

Field 3 - NAME (left adjusted) 

12-character label (columns 12- 
23) for y-axis of multiplot . If 
not specified, the name field is 
selected based on the first 
requested dependent channel. 
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PLOT 


0. EXAMPLE 


1 78 24 88 48 80 71 

1 / / II n i ii j fS I /g j /7 


PLOT 



PLOT 

INDE 2 


PLOT 

SCALE 15 4 


PLOT 

LOG 3 


PLOT 

NAMEEARLY TIMES 

2 

PLOT 

MULT 10 5 

10 . 

PLOT 

CMULT 

4 

PLOT 

CMULT 

21 

PLOT 

CNAMECOINC IDENCE 



10 . 1000 . 

53. 6. 900. 

7 12 6 


Explanation : 

Record Effect 

1 Produce Multiple print and individual plots. 

2 Use Channel 2 for independent variable. 

3 Use plot type 4 (log-log) on the Individual plot of 
Channel 15. Scale the y-axis an the Channel 15 plot 
using 10. as the origin and 1000. as the aaxlRua value. 

4 Plot all nonspecified plots as log-linear. (In this 
case, all except Channel 15). 

5 Label the Channel 2 axes (in this case, the 
independent axis) as "EARLY TIMES.” 

6-8 Produce a aultiple plot of the five dependent channel 

nuabers 4, 7, 12. 6, and 21 versus the independent 

channel nuaber 10. Scale the x-axis using 10. and 53. 
and y-axis using 6. and 900. as the origins and aaxiaua 
values. Note: the plot type will be type 3 as 

specified on Record 4. 


9 Label the y-axis on the aultiplot as "COINCIDENCE." 
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CHAPTER IS - LINEAR SYSTEMS ANALYS IS PROCESSOR ( TRANS ) 


A. FUNCTION 

The Linear Systems Analysis Processor (TRANS) operates on spectral 
density Matrices fro* order 2 up to a aaxiaua of 24. Frequency 
response functions, coherence functions (ordinary, Multiple, and 
partial) and associated confidence liaits are coaputed for single- 
input/single-output and for aul 1 1 pi e- input/s ingle-out put linear 
sys teas . 


B. INPUT/OUTPUT SPECIFICATIONS 
1 . Input 

a . C ontrol Inpu t to Li near Sys teas A nalysis Processor . 

These inc lude : 

o Specification of output channel 

o Plot option for coherence and frequency response 
functions 

o Plot option for confidence liaits 

o Coaputation option for coaplete ordinary coherence 
aatrix 

b . Co ntrol Input s to Data File M a nager (DFM) . FILIN Record 
required to describe the special spectral density foraat 
input file. 

c. Data In put. The special spectral density file as 

generated by either the Time and Frequency Analysis 
Processor (SPEC) or the Power Spectral Density Processor 
( PSD ) . 

I 

2 . Output 

Output from the TRANS Processor includes plots and listings of 
the frequency response (provided in terms of gain in db and phase 
in degrees) and coherence functions. Optionally, the resulting 
transfer functions may be output to SOP. 


C. CONTROL RECORDS 

TRA NS norma ll y re quires very little input, as all the information it 
needs for default execution is on SDF. In summary, the record formats 
are as f ol lows : 
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Record Mo. Plaid 3 


blank 

CONP 

PCON 

OUT 

OROC 

SCAL 

TAPE 


Op era t ion 

Call TRANS 

Set confidence intervals 
Plotting of confidence intervals 
Defines output function 
Ordinary coherence aatrix 
Scaling of plots 

SDP output of the transfer functions 


Foraat 


Record No. 1 ( Optional I 

Default calling of TRANS. Not required If any other TRANS 
record is used. 

Field 1 - TRANS 

Re cord No. 2 (Optional,). 

Used only If it is desired to change the default confidence 
Halts 

Field 1 - TRANS 

Field 3 - CONF , . . 

Field 4 - Confidence interval for coherence functions 

0 ■» 0.95 ( 95% ) is assuaed 

Field 5 - Confidence interval for frequency response 
functions 

0 ■> 0.95 (95*) is assuaed 

Record No. 3 ( Optional! 

This record causes the confidence intervals to be plotted 

Field 1 - TRANS 
Field 2 PCON 

R ecor d No. 4 ( Optional) 

Used only If the output function is other than nuaber one. 
Field 1 - TRANS 

Field 3 - _flUT . , . 

Field 4 - Output channel nuaber (0 ■> function nuaber 1) 
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e . Record No. 5 (Optional) 

Used only if it is desired to print and plot the ordinary 
coherence between all possible input functions in addition 
to the ordinary coherence between the inputs and the 
output . 

Field 1 - TRANS 
Field 3 - ORDC 

f • Rec ord No. 6 (O ption al, up to 25 wav be employed) 

Field 1 - TRANS 
Field 3 - 


Subfield 

Col unns 

Contents 

i 

8-11 

The letters of SCAL . 

ii 

12-15 

GAIN (gain is the only function whose 
plotting way be varied). 

1 i i 

18-17 

One or two digits which indicate 


which data function is to be plotted. 
For exaaple, GAIN 02 would say that 
gain two is to be specially plotted. 
If left blank, all types of functions 
as defined by subfield ii are scaled 
according to this record. 

iv 18-19 Used only if it is desired to define a 

second dependent variable, in which 
case the field contains one or two 
digits indicating the second 
function. 

v 20-23 A single digit to indicate the plot 

type : 

0 . 1 
2 

3 

4 

5 
7 

The balance of the special plotting records is as follows: 


Linear-1 inear 
Linear x, log y 
Log x. linear y 
Log- log 

Linear x, dB y (20 log y) 
Log x. dB y (20 log y) 
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Field 

Co 1 uans 

Linear 

Log 

4 

24-35 

X »1 n 

— 

5 

36-47 

x m 1 1 

nuaber of cycles 

6 

48-59 

Vni n 

— 

7 

60-71 

y»a« 

nuaber of cycles 


g . Record No. 7 ( Opt ionaj j_ 

This is used only if it is desired to output the transfer 
functions onto SDF . 

Field 1 - TRANS 
Pield 3 - TAPE 

If there are aors than twelve input functions, only the 
first twelve transfer functions will be output. Output is 
in the for* of real and iaaginary parts. In orde r for this 
nn t pn t r a be useful, t he "RETAIN" opti on in PSD w_ l 1 1 
probably have to be reque sted if P S D was us ed to se_t_up the 
input aatrix. Otherwise, the nuaber of frequency points 
will not be coaaensurate with MAC/RAN FFT lengths. That 
is, it aay not be possible to inverse transfora the transfer 
functions if the "RETAIN" option is not requested. 

2 . Exaap 1 e 


l 78 24 38 48 80 

I / / i 1 f_3 1 f_4 i fj. i LI i L 

FILIN Filenaae 

TRANS 

END 

Explanation : This is the defaul t execution of TRANS . Filenaae 

is the special spectral density aatrix input to TRANS. The 
standard output will consist of plots and listings of frequency 
response between the output (default is function no. 1) and 
tTie input s in addition to the coherence functions. No SDFs of 
the^output~transf er functions are generated, hence there is no 
need for a FILOUT controller. The END controller terainates the 
run . 
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0. METHOD 


The computational procedures described here will be for computing 
parameters of a mathematical model assuming a p input |x|(t). i * 
1,2,. ...p| and single output |y(t)| linear system. The system 
parameters to be computed are: 


o Frequency response functions between each of the inputs and 
the output . 

o Ordinary coherence functions between the inputs and the 
output, and between all pairs of variables as an option. 

o The multiple coherence function between the output and all 
of the inputs . 

o Partial (conditional) coherence functions between each 
input and the output while conditioning on the other 
i nputs . 


Chapter 9 of DTSA discusses this problem and lists important 
references . 


1 . Comput atio nal Procedures 

The computational procedures necessary for these options 

subdivide into three groups: 

a. Power and cross-spectral density function computational 
rout ines . 

b. A procedure for simultaneously handling p + 1 variables to 
efficiently obtain the spectral density functions among 
all these variables. 

c. The complex variable arithmetic and matrix operations to 
compute the multidimensional linear system parameters. 


The spectral density functions can be computed by the PSD or SPEC 
Processors. In either case a ( p «■ 1 ) by ( p «■ 1 ) set of power and 
cross-spectral density functions is generated and employed as 
the input data for this processor. 


The first operation performed by this processor is a sorting 
procedure. The spectral density functions are each supplied as 
a function of frequency. The Linear Systems Analysis Processor 
eventually must operate on the (p + 1) by (p ♦ 1) spectral density 
matrices, one matrix for each frequency value. 
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# 


The data operated on by the prograa will be a set of spectral 
density matrices at frequencies indexed by k as follows: 


yy (k) 

G y , (k) 

G y2 (k) 

. . . 

G y p ( k ) 

l,< k) 

G t ! (k) 

G i 2 ( k ) 

. . . 

G,p(k) 


G yl ,(k) - 


[ Gp y ( k ) 

Gpi(k) G p 2 (k) 

G p p ( k ) J 


k * 0.1. 

. . . , m 



(15.1) 

The frequency index k 

can represent special 

frequency values 


kf c 

t k - k - 0.1 

m 

■ 

(15.2) 

where f c is the Nyquist cut-off frequency, 
can represent the frequency values 

More generally, k 

f r - f t ♦ kAf k - 0 , 1 

» 

(15.3) 

where 



f t - Beginning frequency 



Af - Frequency increment 


(15.4) 

Note: m might denote the maximum lag value if frequency points 

from (15.2) are used, but might only denote the maximum frequency 
index if (15.3) is used. 

Computation of Frequency Response Functions 




The Mathematical model assumed is indicated in the block diagram 
of Figure 15-1 . 


I The p input variables and the output variable are assumed to be 
f ■ zero mean, stationary. Gaussian processes whenever any 
:1 statistical distribution results are discussed. ‘ rt ie 

'quantities H ly (f), i 1 7 2 p are the frequency response 

function (transfer function) characteristics of the linear 
systems through which the variables are passing to make up y(t) . 
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The variables x,(t) , i - 1,2 p and y(t) are assuaed to be 

discrete (digitized) sequences of N points each. The notation 
for the N discrete points aay be 

x i n - x i ( n&t ) 

X 2 n - X 2 (nAt) 


x p a - x p ( nAt ) 

y „ - y(nAt) n - 0,1,2 N-l (15.5) 

where At is the sampling (digitizing) interval. 

The matrix equation to be solved to determine the frequency 
response function is 


G 1 T 


G 1 i G l 2 ••• G ip 


'1*1, ‘ 

* 

CM 

o 


G z I G 2 z • • • G 2p 


«2, 

• 

m 

. 


• 

Q py . 


. G pl G pl G PP. 


. l*py - 


(15.6) 


where the function argument has been omitted for notational 
simplicity (e.g., Gjj is written instead of G t |(k)). It is 
understood that (15.6) is a function of the frequency index k. 
The matrix and vectors in (15.6) are complex valued and hence 
require complex arithmetic operations for correct manipulation. 
In particular. 

® l i “C|l - JQ t | (15.7) 

where C - C(k) and Q - Q(k) are the appropriate co-spectral and 
quadrature spectral density functions at index value k. 
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The 


solution to (15.6) is 


«lr 


G I 1 G i 2 ... G | p 

- i 

• 

o 

1 

h 2 . 


G l 2 G 2 2 • • • G 2 y 


6 2 y 

, 

* 

. 


. 

.H,y. 


. G pl g p2 G »». 


■ 6 »y. 


(15.8) 


or in simpler notation 

H ly - (G ,,)’ 1 G * y (15.9) 

An individual frequency response function is given by 

l-i 

3 . Computation of Coh erence Functions 

The ordinary coherence functions between the output y and each 
input x, are computed by 

112 2 2 
2 I G I y I G ly +Q, y 

7,y - - i - 1.2 p (15.11) 

G | I G y y G,j G y y 

The multiple coherence function between the output y and all of 
the inputs x t , x 2 x p is computed by 

7y . « * 1 ~ [Gy , G y y ‘ (15.12) 

where G yy denotes the first diagonal element of the inverse 
matrix (G syy )~ 1 associated with G ytI of (15.1). 

Ordinary and multiple coherence functions for the set of inputs 
x t alone are defined by considering the p by p spectral matrix of 
the inputs . 


G 1 1 G ly i - 1.2 p (15.10) 
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G 1 I 
G x z 


G .« 


G j 2 
G 22 


1 p 


G 2 i 


(15.13) 


L “pi 


p2 


' P p 


The ordinary coherence function between any two inputs x, and x ( 
is computed by 


2 

7 i i 


I G i i I 


it “ i i 


2 2 
C j j + Q | i 


ii G i l 


(15.14) 


To obtain the partial coherence function between any input, say 
X|. and the output conditioned on the reaalning (p - 1) inputs. 
G m , is partitioned as indicated below. 


G y x » 


G y y 

G y 1 

1 

G y 2 

Gyp 



mm 

G l y 

G i i 

1 

1 

1 

G 1 2 

... G 1 p 


£ £ 

4— J yy L » y ! 

■ 

G 2 y 

G 2 1 

1 

1 

G 2 2 

G 2 p 

m 

£ £ 

■ 

• 

• 

1 

1 

1 

• 



*— 1 1 ly L — 1 1 1 J 

* 

. G »x 

G P 1 

1 

1 

1 

G p 2 

... G p p J 

- • 

- ' * 


(15.15) 


Then coapute the conditional spectral aatrix 


i, ■ y - y (£ ] 1 £ (15.16) 

<—* yy L — ' y 1 l L — ' llj ly 

This procedure requires the equivalent of the inversion of the 
(p-1) by (p-1) coaplex valued aatrix C t t . An individual eleaent 
G | y | r of the 2x2 aatrix G, y | y can be written in teras of real and 
iaaginary parts as 
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G I y ly * G I y Ip " JQiylp (15.17) 

The partial coherence function between the input x t and the 
output y, conditioned on the other (p - 1) inputs, is now computed 
by 

2 2 

2 I G 1 y Ip I 2 c ly l» + Q|y IP 

7 1 y Ip " - (15.18) 

G l 1 Ip G yy Ip G 1 1 Ip G y y Ip 

Similar results apply for x 2 by interchanging x 2 with x t , for x 3 
by interchanging x 3 with x 2 , etc. 

In the special case of a s ingle- lnput/s ingle-output linear 
system, all coherence functions are identical. This can be 
verified by examining (15.17) for partial coherence and (15.12) 
for multiple coherence. Upon substituting in values when p - 1 , 
these equations will both reduce to (15.11). 

4 . Confidence Limit C ompu tations 

In addition to the computation of the basic parameter estimates, 
the confidence limits for the different coherence functions and 
for the frequency response functions are computed. Confidence 
limits are supplied for gain and phase of the frequency response 
functions . 

The degrees of freedom, n are given by 
n - 28, T 

where B, is the effective spectral resolution bandwidth and T is 
the effective record length. For the special case of the 
s ingle - i npu t / s i ngl e-out put system, p * 1 the equations apply to 
ordinary coherence functions. 

Por a complete definition and interpretation of the confidence 
1 imi ts , see DTSA . 

5 . Number of Digits of P recision Plo t 

Equation (15.18) is not solved explicitly. Rather, a 
simultaneous equation solution is employed. The computations 
are performed in a manner that also produces an estimate of the 
accuracy of results from a numerical viewpoint. In particular, 
it estimates the number of decimal digits which are significant 
in the results. As there is one such result for each frequency 
point, it is essentially a function of frequency for this 
application. A plot is produced of the digits function, and it 
is, in fact, the last plot produced by the processor. 
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